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O Abstract 

(Z2 Active Galactic Nuclei (hereafter AGN) produce powerful outflows which offer excellent conditions for efficient particle accelera- 
C3 tion in internal and external shocks, turbulence, and magnetic reconnection events. The jets as well as particle accelerating regions 
close to the supermassive black holes (hereafter SMBH) at the intersection of plasma inflows and outflows, can produce readily de- 
tectable very high energy gamma-ray emission. As of now, more than 45 AGN including 41 blazars and 4 radiogalaxies have been 
detected by the present ground-based gamma-ray telescopes, which represents more than one third of the cosmic sources detected so 
far in the VHE gamma-ray regime. The future Cherenkov Telescope Array (CTA) should boost the sample of AGN detected in the 
VHE range by about one order of magnitude, shedding new light on AGN population studies, and AGN classification and unifica- 
^s tion schemes. CTA will be a unique tool to scrutinize the extreme high-energy tail of accelerated particles in SMBH environments, 
.' to revisit the central engines and their associated relativistic jets, and to study the particle acceleration and emission mechanisms, 
P^i. particularly exploring the missing link between accretion physics, SMBH magnetospheres and jet formation. Monitoring of distant 
£*-) AGN will be an extremely rewarding observing program which will inform us about the inner workings and evolution of AGN. 
7— I Furthermore these AGN are bright beacons of gamma-rays which will allow us to constrain the extragalactic infrared and optical 
^ backgrounds as well as the intergalactic magnetic field, and will enable tests of quantum gravity and other "exotic" phenomena. 
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1. Introduction 

During the last decade, a new branch of astrophysics has 
emerged at the high-energy end of the electromagnetic spec- 
trum as our cosmos is revealing its amazing richness in the 
TeV gamma-ray band. In the extragalactic domain, more than 
45 AGN have been identified in the very high energy (here- 
after VHE) band above 30 GeV The majority of these sources 
belongs to the blazar class, a peculiar type of AGN with a 
strong Doppler boosted continuum emission emitted by a rel- 
ativistic jet closely aligned with the line of sight. Monitoring 



of bright blazars led to the discovery of flux variability down 
to minute time scales, never so well resolved before in any 
AGN at longer wavelengths, and brought new constraints on 
the physical mechanisms at work |fl] |2) . Another type of ac- 
tive extragalactic sources has recently emerged as TeV emit- 
ters, with the radiogalaxies M 87 H and Cen A (6) for 
which we now have detailed TeV gamma-ray studies. Radio- 
galaxies have a much higher space density than blazars and 
their detection could be less affected by their alignment to the 
line of sight. VHE studies of radiogalaxies are thus likely to 



Preprint submitted to Astroparticle Physics 



April 11, 2013 



contribute significantly to a better understanding of AGN unifi- 
cation schemes. 

The VHE observations already challenge current theories of 
particle acceleration, to explain how particles are accelerated to 
>TeV energies in regions relatively small compared to the fidu- 
cial scale of the black hole event horizon [ 1] [3] [7] [8] • Emission 
models are poorly constrained, with both leptonic and hadronic 
models able to fit most of the available spectral data [9, 10 1 but 
have difficulties to explain fast variability. The energy spectra 
of distant AGN raise specific questions. Their observed shape 
depend on the intrinsic emitted VHE energy spectra, complex 
absorption and/or cascading processes in the AGN host galaxy 
and in extragalactic space, and even on cosmological expan- 
sion and the star formation history of the universe ifTTI . Indirect 
constraints on the infrared diffuse background radiation and the 
stars and galaxies which produce it, can therefore be deduced 
from the observed VHE spectra Ifl2l [T3l [14l . The observations 
taken so far show already that our universe is apparently more 
transparent to VHE gamma-rays than previously thought. Dis- 
entangling the different effects definitively requires high quality 
spectra with better sensitivity and higher spectral resolution at 
energies above TeV energies for large samples of AGN at dif- 
ferent redshifts. 

CTA will provide a unique opportunity to address such is- 
sues. With an order-of-magnitude improved flux sensitivity 
compared to existing instruments, it will offer a large dynamic 
range of more than 10 4 for studies of bright VHE flaring epochs. 
Temporal resolution down to the sub-minute time scales will 
become possible, making CTA a perfect tool for gamma-ray 
timing analysis and for studying AGN micro-variability. About 
four orders in magnitude will be covered in energy, from typi- 
cally 20 GeV to beyond 100 TeV, with an energy resolution of 
typically 10% to 15%, ideally suited for the search for spectral 
features from extragalactic absorption processes. For the better 
defined events with multiple triggered Cherenkov telescopes, 
one arc-minute angular resolution can be achieved, and astro- 
metric positions can be determined with a precision better than 
10 arcsec, which will ensure the reliable identification of AGN, 
and possibly also the detection of VHE emission from compact 
and extended AGN components. 

CTA will be operated as an open observatory, offering a 
multi-functional tool with several configurations and observa- 
tion modes. The flexibility of CTA will be particularly useful 
for the study of AGN. The whole array may be used simulta- 
neously for deep observation of specific fields and targets, or 
some sub-arrays may work independently for monitoring, alert, 
or survey purposes, increasing the global capabilities of the in- 
frastructure. Two CTA sites are foreseen to allow a survey of 
the whole sky, one in the southern hemisphere and one in the 
northern one. The latter will largely be devoted to observe AGN 
and extragalactic sources [15|. 

In this paper we examine how CTA will contribute to an- 
swer open questions on AGN physics such as the nature of the 
black hole magnetosphere, the formation of jets and the accel- 
eration of particles, the total energy budget, and the origins of 
variability, and how it sheds new light on AGN classification 
and unification schemes. After a short presentation of the cur- 



rent knowledge on AGN at VHE in section 2, we analyze the 
performances of CTA in terms of AGN population studies in 
section 3. Section 4 illustrates how the high quality spectra ob- 
tained with CTA can test emission and absorption models and 
constrain various parameters. The capability of CTA to probe 
AGN variability is specifically addressed in section 5. Finally 
we exemplify in section 6 that observing AGN with CTA will 
provide original clues on intergalactic media and diffuse back- 
grounds. 

2. Current status of Active Galactic Nuclei at VHE 

AGN are currently believed to harbor a central massive black 
hole surrounded by an accretion disk of matter spiraling to- 
wards the black hole. In about 10% of the cases (50% for the 
most energetic sources), energetic particle beams are emitted 
along the rotation axis of the black hole, giving rise to the so- 
called "radio-loud" AGN with well-collimated radio jets. There 
is a large variety of AGN classes, with different observational 
characteristics lfT6lfT7l . Standard AGN unification schemes at- 
tempt to classify the sources according to their viewing angle 8. 
The blazar family is the least numerous class, as the blazars are 
believed to be radio loud AGN with their jets aligned to within 
a few degrees of the line of sight to the observer. This leads to a 
strong relativistic boosting effect which amplifies the observed 
luminosity by a factor of ~ £ 4 and shortens the observed flare 
time scales by a factor of 6, where the Doppler factor of the 
relativistic bulk motion is 
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where ft - Vbuik/c, and the bulk relativistic Lorentz factor 
Tbuik - 1/(1 - P 2 ) 1 ^ 2 - This relativistic beaming also contributes 
to the reduction of internal absorption of y-rays by reducing 
the intrinsic luminosity of the source and allowing intrinsically 
larger emitting zones. Indeed, most of the VHE bright AGN 
belong to the blazar family (FigfTJl. Their TeV radiation comes 
from extremely relativistic jets, and strong Doppler boosting 
favors their detection. Blazars include various types of AGN 
such as HBL (high-frequency peaked BL Lac), IBL and LBL 
(Intermediate and low-frequency peaked BL Lac), and FSRQ 
(Flat Spectrum Radio Quasars). Another class of AGN, the ra- 
dio galaxies, has also been detected in the TeV range at low 
redshifts. It is not yet clear how gamma-ray properties of these 
sources, believed to have relativistic jets oriented at larger view- 
ing angles, compare to the ones of TeV blazars. 

The sample of AGN detected at VHE currently includes 49 
published sources, spread in redshifts from z = 0.0018 to z = 
0.536, with 45 blazars (33 HBL, 4 IBL, 4 LBL, 3 FSRQ, and 
one blazar of unknown type), four Fanaroff-Riley type 1 radio 
galaxies, and a few other cases to be confirmed. The additional 
case of the Galactic Center, considered as a weak AGN, will be 
further discussed in section 3.4. Sources can show quiescent, 
low, and highly active VHE flaring states. The AGN fluxes at 
TeV energies range from 0.003 to ~ 20 Crab units, from low 
states to the brightest events. The VHE energy spectra of most 
sources can be described by simple power-laws, with observed 
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Figure 1 : Distribution in galactic coordinates of the VHE active galactic nu- 
clei currently detected by the present generation of Cherenkov experiments 
H.E.S.S., MAGIC and VERITAS, as listed in the TeVCat catalog in October 
2012). The detection of some fifty sources in the TeV domain is now firmly 
confirmed. Various types of AGN are shown, as well as their distribution in 
redshift. 



photon index r o b s between 1.9 and 4.6. Variability is frequently 
observed, namely in about 20 sources out of the 45, despite 
sparse time coverage for many sources. Flux variability has 
been detected on all time scales from years, months, and days, 
down to the minute scale for three flaring sources (Figpl). 

Highly variable VHE events of AGN likely originate at small 
distances (less than 1 pc) from the central engine since rapid 
variability (f var ) on time scales from a few days [ 1 8 1 down to a 
few minutes [40| limits the comoving size of the emitting region 
to 



R' < 



ct var S 



(2) 



where 5, the Doppler factor, is in the range from a few up to 
a few tens, gives R' ~ 10 14 to ~ 10 17 cm. VLBI radio mon- 
itoring of the closest objects confirms this argument by show- 
ing growing evidence of correlation between VHE activity and 
VLBI core evolution EU1 I2T1 l22l l23l . The existence of addi- 
tional extended VHE emission remains an interesting possibil- 
ity mnnanu. 

The number of detected TeV sources per class of AGN ap- 
pears extremely peculiar since the blazars are usually the small- 
est population among all types of AGN seen across the elec- 
tromagnetic spectrum, yet they by far dominate the TeV sam- 
ple. The current TeV HBL sample is not flux-limited and is 
highly biased. The dynamical range of present-day IACT re- 
mains usually below 5000, while active states can easily am- 
plify the fluxes by factors of 10 to 200, and Doppler boost- 
ing by factors of 10 4 to 10 6 , or even more. Comparatively, a 
factor 7 in redshift (from Mrk 421 to 1ES 1011+496 for in- 
stance) would decrease the flux typically by a factor 50 only. 
So the VHE sample appears to be largely incomplete and suffer 
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Figure 2: The big flare of PKS 2155-304 observed in 2006 by H.E.S.S., with 
variability down to a few minute scale |2|. Fig |16| of part 5 shows how CTA 
could have seen it. 



from strong observational biases towards large Doppler boost- 
ing and active states, due to the present sensitivity limits and the 
strategy of observations, often done under VHE and multi wave- 
length alerts. Luminosity functions at VHE are still impossible 
to derive due to poor statistics and to insufficient time coverage 
to firmly distinguish the quiescent stationary states from bright 
flares. 

The spectral energy distribution of TeV AGN appears 
double-peaked, with a first bump in X-rays and a second one 
in gamma-rays (Fig(3]l. Additional emission at lower energies 
comes mainly from the stellar population of the host galaxy 
in the optical range, and from extended jets, hot spots and 
lobes in the radio range. Both leptonic and hadronic scenar- 
ios are widely invoked to describe the radiation processes in 
the VHE domain. The low energy bump that is visible in the 
SED of blazars in the optical to X-ray range is dominated in 
both leptonic and hadronic models by the synchrotron emission 
from a relativistic electron population. In leptonic models one 
generally assumes synchrotron self-Compton processes to ex- 
plain the high energy bump in high-frequency peaked BL Lacs 
(HBLs), with additional contributions from inverse Compton 
scattering on external photon fields to account for the SEDs 
of BL Lacs peaking at lower energies (IBLs, LBLs) and flat- 
spectrum radio quasars (FSRQs). The inverse-Compton inter- 
action (<? + yo — * e + y) °f the X-ray synchrotron emitting 
electrons with ambient photons produces the gamma-ray bump 
at energies hv ~ min[ylhvQ,y e m e c 2 ], where y e is the electron 
Lorentz factor making the peak synchrotron photons. This sce- 
nario is called "synchrotron-self-Compton" (SSC) if the ambi- 
ent photons are the synchrotron emission of the same electron 
population. If the photons come from another radiation field, 
namely from the accretion disk, clouds, or dusty torus of the 
AGN, or from other parts of the jet (e.g. a slow envelope, or up- 
stream or downstream plasma) the scenario is called "external 
inverse-Compton" (EIC). In hadronic models, very high energy 
protons are at the origin of the gamma-ray emission, either di- 
rectly by their synchrotron (or curvature) emission, or through 
their interaction with local gas and radiation background, cre- 
ation of pions and subsequent decay into VHE photons (mainly 
7t° — > 27) with typical energies E y ~ E„/2 ~ 10%E V . De- 
cay of pions into muons also produces neutrinos, and secondary 
electrons which can radiate in the X-ray range. A detailed pre- 



sentation of the above scenarios can be found for instance in 
(271|28l. 
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Figure 3: Example of the double-peaked SED of TeV AGN as seen with current 
instruments, reproduced by a basic SSC scenario. The figure shows the quies- 
cent high energy state of PKS 2155-304 observed during a multiwavelength 
campaign in 2008 with simultaneous data from ATOM, RXTE, Swift, Fermi 
and H.E.S.S. (30). 



An additional effect modifies the shape of the spectra at high 
energy gamma-rays, due to the extragalactic background light 
(EBL) absorption which becomes quite important for high red- 
shift sources. While traveling through the extragalactic space, 
VHE gamma-rays interact with the EBL photons and suffer 
strong absorption due to pair creation, 



TeblTvhe -> e e 



(3) 



with /ivebl x hvvim > (m e c 2 ) 2 . The absorption of VHE pho- 
tons is mainly due to the infrared background generated by the 
integrated light of stars, galaxies and dust (see the EBL article 
in this issue for more details). The observed photon spectrum 
<J>ob s (£\ z) of a source at redshift z is attenuated compared to the 
intrinsic spectrum <$> em (E') emitted at the source, so 



<£ obs (£,z) = e - T '' <£ - ) <D e 



(4) 



where t 7 (E, z) is the optical depth computed along the line of 
sight, and E' = E(l + z). In the energy range 0.2 to 2 TeV, 
the emitted and observed photon indexes T are approximately 
related by 



r obs (z) ~ r em + at(z) 



(5) 



where AT(z) is an increasing function of z |29|. Applying var- 
ious reasonable models of the EBL, one can deduce intrinsic 
photon indices in the range ~ 1.7 to ~ 3 from the observed 
spectra, compatible with data at low redshifts. Future data espe- 
cially from interplanetary space missions should be able to sig- 
nificantly improve the direct measurements of the EBL which 
will be mandatory for a better access to the intrinsic spectra of 
AGN. 

If the intrinsic AGN spectra were all similar, then the distri- 
bution of observed photon index versus redshift should show an 



apparent increase of |r o t, s | with redshift [31], which is not really 
obvious at the moment, even when considering only one spe- 
cific sub-class of blazars (FigHl. This is currently a matter of 
debate, although the departure from standard power laws and 
potential observational biases, such as the current difficulty to 
measure large values of r o t, s , need to be further investigated. 
Non-standard explanations have been proposed, including vari- 
ous evolutionary effects in the blazar population 02113311341 . 
ultra-high energy cosmic rays (UHECR) effects ll35l l36l or 
the existence of axion-like particles which could significantly 
reduce the EBL dimming 11371 [38l . More conservative solu- 
tions exist, such as a hardening of intrinsic spectra during ac- 
tive states 1 39 1, together with the fact that high redshift objects 
might be preferentially detected during flaring states (which is 
indeed the case of 3C 279). However none of them is fully sat- 
isfactory yet. Increasing the statistics on AGN at all redshifts 
and gathering high quality spectra with CTA should clarify this 
intriguing situation. 




Figure 4: The present r„b s -z distribution. This figure includes only published 
sources with confirmed redshift (plus two lower limits at z ~ 0.32 and 0.5) and 
good spectral index determination. The error bars correspond to the statisti- 
cal uncertainties. The observed spectral indices are mostly derived in the 0.2-2 
TeV range, the relevant EBL range in that case being 0.25-2.5 eV. The spectral 
energy distribution of the EBL departs from a power-law in this range, show- 
ing an emission bump. Bracketing the spectral photon number density versus 
energy by power laws, De Angelis et al (2009) |37j derived an approximate 
analytic expression for the optical depth within the Franceschini et al (2008) 
model |41 1 in the local Universe. The two grey curves in the figure correspond 
to the expected lower and upper limit after taking into account the cosmological 
effects, assuming that T ~ 2.35 at z = 0. Apart from two radiogalaxies at very 
low redshift, various types of blazars are shown in color (red = HBL, black = 
LBL, blue = IBL, green = FSRQ). Triangles indicate confirmed flaring states. 
A few sources are shown in both flaring and non-flaring states (Mrk 501, PKS 
2155 and 1ES2344), and two are shown only in flaring states (W Comae and 
3C 279). 



3. AGN population studies with CTA 



The jump of sensitivity with CTA will offer large samples 
of VHE sources and open the way towards statistical studies 
of the blazar and AGN populations. Large and homogeneous 
samples of sources are necessary for unbiased statistical ana- 
lyzes and to construct useful luminosity functions. To be able 
to gather them, a significant step forwards will be to firmly de- 
tect and study quiescent stationary VHE states for comparison 
with bright flares of blazars and AGN. Currently, this appears 
possible only for a few bright sources such as the BL Lac PKS 
2155-304. Dealing with real stationary fluxes (if any) will at 
last allow a good statistical approach to the gamma-ray AGN 
samples, by comparing objects in the same activity state. This 
will become possible with CTA. It will allow us to explore the 
relation between HBL, IBL, LBL and FSRQ, and to clarify the 
validity of the so-called blazar sequence and its extended ver- 
sions 22] 03] 04] |45| as well as the recent "blazar envelope" 
view 11461 or other tentative unifying schemes which try to ex- 
plain the trends by leptonic scenarios with a decreasing rela- 
tive importance of the SSC to the EC emission from HBL to 
FSRQ, or to reproduce all observed SED as a function of black 
hole mass and accretion rate Il47ll48l . Furthermore, comparing 
properties such as variability characteristics between quiescent 
and flaring states will contribute to elaborate a global scenario 
for VHE phenomena in AGN. Is it the same emission mecha- 
nism which dominates both quiescent and active states, or are 
there quite different origins for the VHE emission, with sev- 
eral mechanisms at work depending on the activity level of the 
source? CTA will answer such questions, which is important for 
a better understanding of the AGN physics and also to elucidate 
the observational biases known to affect their VHE detection. 

The distribution in redshift of known BL Lac sources peaks 
around z ~ 0.3, with a large majority of the population within 
Z < 1 l49ll50l . a domain of redshift easily reachable with CTA 
in flaring and even non-flaring states for the brightest sources 
(see Fig|5]and FiglS). Conversely, bright FSRQ should be de- 
tectable beyond redshift z ~ 2 with long exposure time of 50 
hours (see Figj6] and Figl9]». CTA will therefore clarify the 
actual distribution of observed photon index versus redshift 
(FigJ4j) and should offer the possibility to study evolutionary 
effects at VHE, at least for the blazar class. 

3.1. The population of blazars 

What will be the best strategy to gather complete samples of 
sources? A blind imaging survey is the most unbiased way to 
observe in a given waveband, and potentially includes the possi- 
bility of discovering unexpected objects, such as isolated black 
holes expelled from host star clusters IBTl |52l or final evapo- 
ration of primordial black holes l53l [54). Given our present 
ignorance of AGN properties and statistics at VHE, it is not yet 
certain how many AGN can be detected by deep surveys for a 
realistic observing time with CTA. 

Estimates have been done for the blazar class, including BL 
Lac sources and FSRQs, based on a blazar gamma-ray lu- 
minosity function model consistent with recent Fermi results 
[55, 56 1 and under the conservative assumption of the standard 
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Figure 5: Significance of detection in units of the standard deviation a versus 
redshift for HBL flares of the type of PKS 2155-304, assuming a typical flare 
with an integrated flux of 10 Crab atz = 0. 1, a spectral index of 2.5, and 3 hours 
of observing time, for three different CTA arrays. 
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Figure 6: CTA detection at highest redshifts: the Fermi blazar 
2FGLJ1504.3+1029 at redshift z = 1.839 as possibly seen by CTA. Here the 
VHE spectrum of the blazar during a quiescent state has been extrapolated from 
the Fermi data assuming a log-parabolic spectral shape, and including EBL ab- 
sorption effects deduced from [41], which are clearly visible in this figure. De- 
tection by CTA can be achieved in 50 hours at 13.6<x (array configuration B, 
zenith angle of 20 degrees). A small sample of blazars should be reachable by 
CTA at redshifts ~ 2 during quiescent states. Bright and long flaring states may 
allow to obtain useful spectra from sources at even higher redshifts, such as the 
distant Fermi blazar 1FGL J1344.2-1723, at redshift z ~ 2.49. 



blazar sequence to simulate the SED. The number of blazars 
detectable by CTA in blank fields has been estimated for var- 
ious array configuration (see also the survey article in this is- 
sue). The best case leads to about 0.36 blazars detected above 
30 GeV in 40 deg 2 , for 50 hours of exposure time. This results 
into 370 blazars potentially detectable by CTA, for an all-sky 
survey with 50 hours/FoV. This requires a very long observing 
time, and more than 30 years [J comparable with the expected 



Throughout this section, we assume a CTA field of view (FoV) of 7 degrees 
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Figure 7: Cumulative redshift distribution of blazars above 30 GeV for various 
CTA array configurations (50 hours observing time for each field, 20 degrees 
zenith angle observations). Here SED have been simulated following the stan- 
dard blazar sequence scenario, which seems to underestimate the importance of 
VHE emission, at least at high redshifts (see Fig|6]and Fig|9]for comparison), 
so these distributions should be considered as lower limits. 



lifetime of CTA. However, relatively bright objects should read- 
ily come out in a full-sky survey at ~ 1 % Crab level in about 
one year. The expected cumulative redshift distribution in the 
entire sky is shown in Figj7] For a conservative view of blazar 
statistics, which tends to minimize VHE emission, CTA has 
the potential to find at least 20 blazars at redshift larger than 
1 with its configuration B [57. 58|. Should the standard blazar 
sequence be revised with more high redshift TeV blazars, still 
higher detection rates would be expected. 

Blank field sky surveys are mandatory because they provide 
samples of sources with minimal observational biases, espe- 
cially from other waveband observations or from detections 
triggered by alerts, which is essential for any reliable statisti- 
cal analysis of the blazar population. As deep surveys are very 
time consuming, they will be extended over many years, in the 
fields of specific targets. One can anticipate that for a FoV of 
7 degrees, more than 100 serendipitous discoveries of blazars 
should be obtained in the field of other extragalactic sources in 
about ten years. 

3.2. Extrapolating from Fermi sources 

Follow-up observations of sources already detected at lower 
frequencies offer another promising approach — biased but 
faster — in the search for TeV AGN. Here we consider the case 
of the Fermi satellite, which has already provided detections of 
a large sample of AGN in the high energy range. The detection 
of blazars in the Fermi LAT energy band (~0.1 to 100 GeV), 
i.e. at energies below or partially overlapping with the ones 
accessible with CTA, provides important information about the 
intrinsic spectral properties of extragalactic sources, which are 
mostly unaffected by y-ray absorption on the EBL in the Fermi 



and 1500 hours of observing time per year. One can anticipate that typically 
one third of the whole CTA observing time, including the two sites, could be 
devoted to AGN programs 



LAT band. The combination of same-epoch Fermi and CTA 
measurements of AGN spectra will be of great importance since 
it will provide a broad-band measurement over ~ 5 decades in 
energy. 

In order to make predictions for CTA, we extrapolated 
Fermi/LAT AGN spectra into VHE frequencies taking EBL at- 
tenuation into account [41 1. Two independent analyses, based 
on slightly different selections of sources and assumptions on 
their detection, lead to compatible estimates and predict a large 
sample of AGN within the reach of CTA. 

For this work, the effective areas and cosmic ray backgrounds 
simulated for different array configurations have been used to 
estimate the instrumental response of CTA (cf. the Monte Carlo 
contribution in this issue). For simplicity, we consider here full 
CTA arrays in north and south, although a full configuration 
with arrays like E, I and C is not planned at the moment for the 
northern site. With 24 months of accumulated data, the 2FGL 
catalog contains 1873 point sources characterized in the 100 
MeV to 100 GeV energy range [59|. Many of these sources 
have not yet been discovered in VHE due to the lack of sensitiv- 
ity of existing instruments, but they will be accessible to CTA. 
In fact, about 85% of the VHE active galactic nuclei detected by 
ground-based Cherenkov observatories are found in the 2FGL 
catalog and an extrapolation of the 2FGL data to higher ener- 
gies seems a sensible procedure for the compilation of a mock 
catalog of CTA sources. 

All associated/affiliated extragalactic sources were selected 
from the 2FGL catalog [59. 60]. Out of these 1098 extragalactic 
sources, only those with a measured redshift can be used for 
our studies, as it is needed for applying the EBL absorption. 
For a first study, aimed at providing a conservative estimate of 
source counts, AGN were only selected for further processing 
if no analysis problems were flagged in the 2FGL catalog. In a 
second study, designed to include a maximum number of Fermi 
sources, no cut was applied to the 2FGL flags. Technical details 
on the source selection and the significance estimates can be 
found in |Appendix A| 

The first study focuses on the array configuration B with the 
best sensitivity at low energy. It yields about 170 detectable 
AGN within 50 hours of maximum exposure time per source 
(Fig. [8] panel b). Out of these, 27 sources (resp. 70) should 
already be detectable after half an hour (resp. 5 hours) of expo- 
sure time with the full array (Fig. [8] panel a). With a maximum 
exposure time of 150 hours per source, which is not unusual for 
observations of very faint sources with the current IACT, about 
230 sources should be detectable (Fig. [8] panel c) in less than 
10 years. 

Apart from augmenting the total number of extragalactic TeV 
sources, CTA should also increase dramatically the number of 
objects visible at high redshifts. Fig. [9] shows the distributions 
of AGN as a function of redshift. The most distant quiescent 
AGN predicted here is at z — 2.2 and under certain conditions 
this limit might be raised to even higher redshifts. In particu- 
lar, our estimations indicate that certain flaring FSRQs with a 
gamma-ray flux increase of a factor of 10 and moderate spectral 
hardening AF = 0.3 could produce further detections at z > 2. 

The second study explores the performances of different ar- 
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Figure 8: Skymaps of AGN with confirmed redshift from the second Fermi 
catalog (2FGL) that are potentially detectable by CTA in 5 hours (panel a), 50 
hours (panel b) and 150 hours (panel c) maximum exposure time per FoV (as- 
suming the array configuration B and a 20 degree zenith angle over the whole 
sky). It should be noted that several extragalactic sources that are already 
known to be TeV emitters do not appear in this figure, due to the tight se- 
lection criteria applied here and the existence of non-Fermi VHE sources. Such 
skymaps should be obtained in less than two months (a), in about three years 
(b) and in less than 10 years (c) with CTA, assuming 1500 hours observing time 
per year. 



ray configurations (see TablefTll. We obtain typically more than 
> 140 extragalactic detections, for 50 hours of maximum ex- 
posure time per source. In particular, we note that CTA will be 
most efficient for hard sources (T < 2) and shall reveal all the 
complexities of the extragalactic population extending beyond 
the reach of current detections. For softer sources that are cur- 
rently only accessible during flares, CTA will provide unique 
access to quiescent states. 

The array configurations that do not include large-size tele- 
scopes (LSTs), such as array "C", yield the poorest perfor- 
mance. They give access to significantly fewer FSRQs, which 
usually have very steep spectra, but the same is also true for BL 
Lacs. Array B, with the best coverage at low energies, yields 
the best results in terms of source statistics. The "compromise" 
solutions, such as configuration E, and the northern array NA, 
remain a good option. 

The results from these two studies illustrate the remarkable 
capabilities of CTA compared to current IACT The detection 
of at least 140 extragalactic sources is expected in less than 
two years with a maximum of 50 hours per source. Actually, 
the total number of Fermi AGN detectable in 50 hours would 
reach about 370 if the 2FGL BL Lacs without known redshifts 
were detected with the same proportion (namely one third for 



Figure 9: Redshift distributions of the AGN with confirmed redshift from the 
second Fermi catalog (2FGL) that are potentially detectable by CTA in 5 hours 
(panel a), 50 hours (panel b) and 150 hours (panel c) maximum exposure time 
per FoV (assuming the array configuration B and a 20 degree zenith angle over 
the whole sky). These distributions correspond to the skymaps shown in Fig|8] 
The same colour code is used: brown for BL Lacs, green for FSRQs, red for 
radio galaxies and magenta for other types of AGN. 



our two studies) than the ones with confirmed redshift q] In 
fact, the main difficulty in elucidating the full BL Lac popula- 
tion lies in obtaining direct redshift measurements from their 
mostly featureless optical/UV spectra. A possible workaround 
might come from a direct measurement of the shape of the EBL, 
which would allow us to set an upper limit on unknown red- 
shifts (see EBL contribution in this issue). The artificial break 
introduced here at 100 GeV for hard sources might also lead to 
an underestimation of the number of detections, but on the other 
hand, some of the sources might have intrinsic spectral breaks 
or cutoffs above the Fermi LAT energy range, which would re- 
duce their signal in the VHE band. Only CTA will inform us on 
the actual spectra beyond 1 TeV for most of these sources. 

Discoveries such as extreme AGN or extragalactic sources 
not adapted to the Fermi LAT band and with emission peaking 
in the CTA energy range should further increase the sample of 
AGN seen by CTA. Given that currently 6 out of 45 TeV AGN 
have not yet been detected with the Fermi LAT [62], one can an- 
ticipate a fraction of about 15% of missing CTA sources when 
making predictions based on the 2FGL catalog. All these es- 
timates are in agreement with the lower limit predictions made 
from the blazar population analysis presented in part 3.1 (which 



However, the percentage of detection could be smaller for BL Lacs without 
redshift, since they have a higher probability to be remote objects 



Table 1: Number of detectable Fermi AGN with redshift for different array configurations (50 hours of maximum exposure time). AGN with unknown type are 
classified as "other AGN". 
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includes Fermi and non-Fermi blazars, with or without known 
redshift). They prove that within its lifetime, CTA should pro- 
vide significant samples of AGN of various types, suitable for 
statistical analyses of the VHE population. 

3.3. Radiogalaxies and extended sources 

Detection and monitoring of four radio-galaxies at TeV ener- 
gies definitely proved that AGN other than blazars actually ra- 
diate at VHE. This is a key information which came out rather 
recently. Indeed, although quite rewarding by itself, detecting 
only blazars at VHE energies provides biased information about 
AGN jets, always seen at small viewing angles. Studying AGN 
with moderate or negligible Doppler boosting with CTA will 
provide access to a 2D view of AGN jets at VHE and will shed 
new light on our current understanding of radio-loud AGN. In 
particular, the performances of CTA should allow a detailed 
analysis of the relation between VHE and non-thermal radio 
emission and will contribute to the long-standing question of 
the origin of radio-loudness in AGN. 

The four radiogalaxies detected so far at VHE, namely M 87 
0|551, Cen A |5J, IC 310 |5J|55j, and NGC 1275 |>51, have 
all been tentatively classified in the literature as Fanaroff-Riley 
type 1 radiosources, each with some "peculiarities". Apart 
from that, these active galaxies and their nucleus and jets show 
very different properties which, based on knowledge already 
gathered at lower energies, do not suggest any specific promi- 
nent common features, except being TeV sources. These VHE 
galaxies form an emerging class of AGN, which should find 
its place within AGN grand unification. This simple fact illus- 
trates that observing at VHE probes new aspects of AGN not 
yet explored at any other energies, and offers a fully indepen- 
dent tool of investigation. One possibility would be that the 
VHE band directly catches the emission from the base of an in- 
ner beam or jet when it has a right orientation, independently 
of any other properties and classification of the radiosource at 
others frequencies. 

This emphasizes the strong interest of studying such types of 
sources with CTA, but precludes at the moment any convincing 
prediction of a sample to be detected at VHE in the future. In 
particular, in the current small sample of VHE radiogalaxies, 
the average detected TeV fluxes do not seem to be related to 
the average non-thermal radio and X-ray fluxes of the sources. 
One common trend is that the four sources are located in rich 
environments, and show sign of galaxy interaction or mergers. 



Both M 87 and NGC 1275 are dominant cluster galaxies with 
a very massive central black hole. Cen A is presumably a re- 
cent merger, in a group of galaxies, and IC 310 is located in 
the Perseus cluster where radio jets can interact with the in- 
tracluster gas. Other common properties could be to have an 
intermediate viewing angle (about 20 degrees for M 87, 40 for 
Cen A, between 20 and 50 for NGC 1275 and < 38 for IC 310), 
weak or moderate Doppler beaming, and no direct alignment 
between their radio compact VLBI core and their extended ra- 
dio structures (grossly misaligned by ~ 70 degrees in M 87 and 
~ 45 in Cen A), possibly enhanced by projection effects and in- 
homogeneous external medium, or related to specific properties 
of their central engines. Moreover, the detection of some tran- 
sient BL Lac-type phenomena has been reported or discussed 
in the literature for the four radiogalaxies discovered up to now 
at VHE. Such properties should help to identify promising can- 
didates for further observation with CTA. 

Among this new class of AGN, M 87, the first non-blazar 
detected in the VHE range, has been the most studied in the 
literature. It questions our global understanding of TeV emis- 
sion scenarios for AGN, and to some extent our general view of 
AGN classification. Even assuming a high Doppler factor, the 
TeV variability of M 87 requires very small emitting zones, of 
the order of a few Schwarzschild radii R s of its 3 x 10 9 M Q black 
hole, under causality argument 01141151 12111 . This raises the crit- 
ical question of particle acceleration mechanisms in such small 
regions, and excludes the Virgo cluster, the radio lobes, the host 
galaxy, the large scale jet and its brightest knot A as dominant 
TeV emission zones. Three different main emitting zones have 
been considered for M 87, (i) the peculiar knot HST-1 located at 
about 65 pc from the nucleus, (ii) the inner VLBI jet, and (iii) 
the central core itself, namely the accretion disk or the inner 
black hole magnetosphere. 

The coordinated campaign in 2005 found a possible corre- 
lation between a VHE flare and an X-ray outburst of HST- 
1. However, another multi wavelength campaign organized in 
2008 concluded that the X-ray light curve of HST-1 obtained 
by Chandra does not follow the VHE one. Conversely, the ra- 
dio and X-ray emissions from the core are correlated with the 
VHE flux. In radio, regular monitoring of M 87 by the VLBA 
at 43 GHz also allows one to explore the sub-mas scale in order 
to probe the jet formation and collimation zone at about 100R S 
from the black hole, with an outstanding angular resolution of 
0.21 mas x 0.43 mas (0.5mas ~ 0.04pc ~ 140R S at the distance 



of M 87). The VLBA detected a significant rise of the flux of the 
radio core at the time of the 2008 VHE activity, together with 
enhanced emission along the VLBI jet. These results favour 
scenarios where most of the VHE emission comes from the in- 
ner VLBI jet (multi-zone models inspired from standard blazar 
scenarios [67, 68 69|) or from the central core (particle accel- 
eration in the black hole magnetosphere 1170117111721 ). However 
the situation still remains unclear, and suggests the existence 
of different types of VHE flares, as discussed in a recent work 
reporting on the 2010 large joint monitoring campaign [73 1 . 

Current magnetospheric models for M 87 ll70l ITTI l72ll imply 
that minimum variability timescale should always be larger than 
a few Schwarzschild light crossing times R s /c ~ (0.2-0.4) days 
(in the absence of strong Doppler effect) and that the TeV spec- 
tra should exhibit a clear break, due to internal yy absorption 
and maximum energy constraints, well below 50 TeV. Confirm- 
ing or rejecting such magnetospheric scenarios with "particle 
acceleration close to the SMBH" by future high-sensitivity ob- 
servations with CTA would boost our understanding of AGN 
central engines: in the former case, new jet physics beyond 
current developments would be demanded, in the latter case, 
a strong link between jet formation, particle acceleration and 
disk physics would be established. 

In the second radiogalaxy discovered at VHE, the nearby 
source Cen A, the origin of the dominant VHE signal is even 
less clear than in M 87, as both the radio core and the kpc 
jets of Cen A are within error bars of its position on the sky. 
This results in many possible emitting zones such as the black 
hole magnetosphere, the base of the jet, the large scale jets 
and inner lobes, or even a pair halo in the host galaxy. A bet- 
ter accuracy on the absolute astrometry at VHE expected with 
CTA and obtaining high quality light curves with good tempo- 
ral coverage should clarify this decisive issue, and at least dis- 
tinguish between a dominant core emission or a dominant ex- 
tended component related to the kpc jets. Indeed, very extended 
diffuse gamma-ray emission has been recently found by Fermi 
at lower energies [74- 1 (see FigfTO}. The gamma-ray emission 
above 100 MeV coming from the giant lobes can be described 
by EIC models on the cosmic microwave background and the 
EBL. This emission is quite important, with a total flux slightly 
higher than the one from the core [75] and a power comparable 
to the kinetic power required in the jets. The high sensitivity 
and improved angular resolution of CTA should allow one to 
look for and possibly map any extended structure at higher en- 
ergies, thereby opening a completely new view on VHE particle 
acceleration, transfer and radiative losses. A simple extrapola- 
tion of the Fermi halo to the VHE band shows that its VHE 
counterpart remains out of reach by CTA, assuming the same 
spatial extension of about 2 degrees. However, VHE detection 
can be expected from regions of enhanced gamma-ray emission 
in shocks, knots, or hot spots, which could be identified with the 
angular resolution of CTA. Conversely, the detection of an ex- 
tended VHE halo around M 87 would be very challenging but 



could be possible with CTA (see Fig 111, assuming, as seen by 
Fermi in Cen A, a flux in the lobes of 50% the flux detected in 
the TeV range during a low state of the source [ 3 ] for an ex- 
tension of 0.2 degrees corresponding to the extension of radio 



maps. 

In contrast to the three other VHE radiogalaxies, IC 310 was 
initially not recognized as having a specifically remarkable non- 
thermal activity. Its serendipitous discovery at VHE in the field 
of NGC 1275 emphasized our current poor knowledge on VHE 
populations. Only recently it appeared in the Fermi catalog and 
was then identified as a potential TeV source. Two emission 
zones can be considered, the central engine and the inner jet 
as commonly described for TeV BL Lacs, or the bow shock 
created by interaction of the fast moving host galaxy with the 
intracluster gas l65l . Astrometric and angular resolution ca- 
pabilities of CTA should distinguish between them. However, 
the first option appears favored because of the detection of few- 
dayscale variability. Indeed, IC 310 was already mentioned in 
the literature as a FR I source which may have a non-thermal 
activity related to the BL Lac phenomena, but at weaker levels 
than characterized by the standard definition of BL Lac objects 
[76]. Moreover recent VLBI data show a blazar- like one-sided 
core-jet structure at intermediate angle to the line of sight [771 ■ 
VHE instruments are therefore possibly on the way to solve the 
long standing problem of the "missing BL Lac" and to firmly 
identify the still elusive transition population between beamed 
BL Lacs and unbeamed FR I galaxies [78], a difficulty of the 
standard unification scheme which proposes that BL Lac are 
FR I radiogalaxies seen along their jet axis. Surveys at VHE 
could have the capability to recognize a population of low lu- 
minosity or misdirected BL Lacs, difficult to identify at lower 
energies, and thus "bridge the gap" between genuine BL Lacs 
and FRI radiogalaxies. It will be interesting to further investi- 
gate such a view in the context of recent blazar classification 
scenarios [50|. 

Generally speaking, nearby radio galaxies offer the oppor- 
tunity of unique studies of extreme acceleration processes in 
relativistic jets and in the vicinity of supermassive black holes. 
Given the proximity of the sources and the larger jet angle to 
the line of sight compared to BL Lac objects, the outer and 
inner kpc jet structures are potentially resolvable by CTA, en- 
abling us to look for possible VHE radiation from large scale 
jets and hot spots besides the central core and VLBI jet, and 
to spatially pin down the main site of the emission. With the 
help of simultaneous multiwavelength observations and tempo- 
ral correlation studies, different sections of the jet and the core 
can be probed, down to the smallest pc (milliarcsecond) scale, 
only accessible to VLBI radio observations or timing analysis. 
Further studies of variability with CTA will strengthen the lim- 
its on the size of the emission region and clarify the correla- 
tions with other wavelengths. Long-term monitoring and the 
search for intra-night variability would be two major goals to 
constrain the physics and start characterizing this new popula- 
tion of sources. Remembering the basic classification of extra- 
galactic radio sources, one could consider highly variable VHE 
radiogalaxies as likely core-dominated gamma-ray sources, and 
poorly variable ones as possibly lobe-dominated gamma-ray 
sources. This VHE population is still lacking a standard uni- 
fying model and deserves further analysis. 

Recent and fast developments on VHE radiogalaxies show 
that present VHE instruments start to provide an original view 
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Figure 1 0: The diffuse gamma-ray emission detected by Fermi from the giant 
lobes of the radiogalaxy Cen A 1741 . Depending on the actual VHE spatial 
distribution, and the sensitivity and angular resolution performance of CTA, 
the structure of the extended VHE emission on the kpc scale in the central part 
of the galaxy can be probed (here the white circle corresponds to the LAT PSF 
of 1 degree). 
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Figure 1 1 : CTA detection of extended extragalactic VHE emission: under cer- 
tain conditions, a halo of extended TeV emission around M 87 could be within 
reach of CTA. The figure shows the expected spectrum from the CTA array con- 
figuration I, for 100 hours of observing time, assuming an extended emission 
of 0.2 degrees with a flux of 50% of the total flux detected during low states of 
M 87 (see text). 



of non-thermal activity from central AGN engines and inner 
jets, with the capability to directly probe a very specific re- 
gion, still not fully identified and unreachable by other means, 
in the close vicinity of SMBH, such as for M 87. The next 
generation of IACT will explore this still missing link between 
SMBH magnetospheres and the physics of jets and extended 
radiosources. One can also anticipate that it could provide de- 



cisive constraints on the fundamental question of the total en- 
ergy budget of some non-thermal sources where the contribu- 
tion of the extended gamma-ray emission appears quite signifi- 
cant, such as for Cen A. 

3.4. Seyfert galaxies 

There is a growing evidence that relativistic jets are not only 
seen in blazars and radio galaxies but in several types of Seyfert 
galaxies as well. About 5% of narrow-line Seyfert 1 (NLS1) 
galaxies are radio-loud (RL) f79l . and show fiat spectra to- 
gether with variability in the radio band, suggesting the pres- 
ence of relativistic jets. This hypothesis has recently been 
confirmed by the detection of a small number of RL-NLSls 
with Fermi-LAT 11801 . The measured GeV spectra are typi- 
cally steep, with F = 2.5 - 2.8, which makes the detection of 
RL-NLSls with CTA challenging. However, at least two RL- 
NLSls (PMN J0948+022 and SBS 0846+513) have shown sig- 
nificant variability ll8Tll82l . with gamma-ray luminosities that 
would make them detectable with CTA during gamma-ray out- 
bursts. 

The broadband emission of RL-NLS1 galaxies appears sim- 
ilar to that of blazars, being explained by SSC and EIC scenar- 
ios. Thermal emission from the disk, broad-line region (BLR) 
and the infrared torus dominate the IR/optical output, and pro- 
vide seed photons for EIC scattering |80|. Even though the 
radiative properties of the jets in RL-NLS 1 galaxies look sim- 
ilar to that of blazars, they originate in different host galaxies. 
Blazars are hosted in elliptical galaxies, while RL-NLSls are 
likely in spirals. In addition, gamma-ray emitting RL-NLSls 
have inferred SMBH masses typically 1-2 orders of magnitude 
smaller than blazars, while their accretion rate reaches extreme 
values, up to 80% of the Eddington rate, which have never been 
found in gamma-ray loud AGNs but are usual for NLSls [83 1. 
Detecting RL-NLSls with CTA would significantly increase 
the range of parameters of the accretion-ejection process ex- 
plored at VHE. 

Two Seyfert 2 galaxies have also been detected by Fermi- 
LAT, namely NGC 1068 and NGC 4945 (84). Both galaxies 
exhibit AGN and starburst activity [85, 86 1, and a simple extrap- 
olation of their GeV spectra yields VHE flux levels detectable 
with CTA. Given their similarity with the known gamma-ray 
emitting starburst galaxies, it seems likely that their GeV emis- 
sion originates in the interstellar medium of the galaxy. How- 
ever, 1 84 1 find evidence that the AGN component could be dom- 
inating in NGC 1068. A detection of a significant spectral 
steepening or flux variability in the VHE range would point to- 
wards AGN emission, disentangling it from the steady starburst 
component. 

Classical radio-quiet Seyfert galaxies could also emit 
gamma-rays originating from starburst or AGN activity. Clas- 
sical Seyferts also show jet-like structures, but with flows that 
are typically slow, weak, and poorly collimated when compared 
to the relativistic jets of radio galaxies and blazars. So far, 
only two possible associations between GeV sources and radio- 
quiet Seyfert galaxies have been found (ESO 323-G077 and 
NGC 6814) but chance spatial coincidences with these objects 
cannot be ruled out [ 87 1 . 
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A detailed characterization of the gamma-ray emission from 
the different classes of Seyfert galaxies will test our knowledge 
of jet launching mechanisms in AGNs with SMBH masses and 
accretion rates very different from the well-studied gamma-ray 
loud blazars. This will shed light on fundamental questions in 
AGN physics like formation and propagation of jets, the phys- 
ical cause of radio-loud/radio-quiet distinction, the fundamen- 
tal parameters governing the central engine, and how the host 
galaxy influences the active nucleus. 

3.5. Low-luminosity AGN, supennassive black holes, and the 
Galactic Center 

Applying standard SSC scenarios to other low luminosity 
AGN (LLAGN) such as NGC 4278 shows that they may be 
detected by CTA [88]. Furthermore, scenarios developed for 
M 87 show that magnetospheres of rotating SMBH can easily 
generate VHE particles and radiation, as long as the accretion 
disk remains under-luminous to avoid strong internal absorp- 
tion by the ambient radiation. This raises a very general ques- 
tion. Is VHE emission a generic feature of AGN and SMBH? 
Indeed, many normal galaxies harbor SMBH in their nuclei and 
the search of some VHE signal from them will be a challenging 
topic for CTA which could shed new light on the crucial ques- 
tion of missing SMBH, AGN evolution and feedback between 
AGN and host-galaxies. Spheroidal systems such as elliptical 
or lenticular galaxies, and bulges of early-type spiral galaxies, 
are believed to host SMBH with masses between 10 6 and 10 9 
solar masses. During the early stages of galaxy evolution, the 
SMBH accrete matter at high rates and are observed as bright 
QSOs. The radiative output in the optical decays from redshift 
Z > 3 to z - by almost 2 orders of magnitude. The majority 
of SMBH in the local universe are hosted in those evolved sys- 
tems with low accretion rate, not embedded in dense radiation 
fields. This enables high energy and VHE y-rays, if generated, 
to escape from the nuclear region without suffering from strong 
absorption via photon-photon pair absorption. A good candi- 
date for such detection could be, for instance, NGC 1399, as 
discussed in detail in [89, 90]. Any positive signal of this kind 
would be a real breakthrough for AGN and SMBH physics as 
well as for co-evolution schemes of SMBH and galaxies. It 
could also raise important issues related to the origin of extra- 
galactic cosmic rays, and to the still debated results from the 
Pierre Auger Observatory (PAO) [91 1, and offer an explanation 
for the correlation claimed by the PAO collaboration between 
the arrival directions of the most energetic cosmic rays and the 
spatial distribution of low redshift AGN, which mainly involves 
apparently weak AGN. Further statistical analysis and monitor- 
ing of AGN samples at very high energies will tackle these still 
open questions. 

One remarkable case of weak AGN is the Galactic Center. 
Indeed, TeV emission has been detected from the central region 
of our Galaxy, although it is not yet known whether it comes 
from the SMBH Sgr A* itself or from the pulsar wind neb- 
ula G359.95-0.04 at 8.7 arcseconds, still inside the spatial error 
bars of present IACT, but likely resolvable with CTA. More- 
over, an interesting strategy to fix the position of the VHE emis- 
sion at the milliarcsecond scale is to look for pair-production 



VHE eclipses, when stars orbiting the central SMBH approach 
the line of sight. Such phenomena can be very well described 
since trajectories are precisely known for several stars. Light 
curves can be predicted and appear within reach of CTA if the 
VHE emission zone is really compact |92|. Any positive de- 
tection of such new-type of eclipse in the next decades would 
allow to reconstruct for the first time the position of a VHE 
source with a tremendous milliarcsecond accuracy. 

If CTA concludes that the VHE source in the Galactic Cen- 
ter region is actually associated with the central SMBH Sgr A*, 
one could then try to extrapolate its characteristics to the galac- 
tic nucleus of other nearby galaxies, such as M31 (Andromeda). 
Although its distance (770 kpc) is about 100 times larger than 
that of Sgr A*, its nuclear X-ray luminosity (L x ~ 10 36 erg/s 
[93 1) is about 1000 times higher (L x ~ 10 33 erg/s for Sgr A* 
|94|). Therefore, crudely assuming Lvhe k Lx, one can ex- 
pect a VHE flux from M31 of about 10% of the one of Sgr 
A*, which would be within reach by CTA. This is an additional 
clue for future possibilities of studying weak SMBH in nearby 
galaxies. 

Transient phenomena may also offer new signatures of 
SMBH, like the recent event (Swift J164449.3+573451) dis- 
covered by Swift in X-rays which may reveal a Tidal Disrup- 
tion Event [95 , 96 1. Indeed, the tidal disruption of a star by the 
otherwise "dormant" SMBH of the compact galaxy hosting this 
burst may have activated the formation of a beamed jet. De- 
pending on the synchrotron or Inverse-Compton origin of the 
X-ray emission, high energy gamma-rays could be produced 
1971 . VHE instruments should therefore contribute to constrain- 
ing such phenomena that are still basically unexplored, which 
may bring to light the basic building block of the accretion- 
ejection cycle around SMBH. 



4. Particle acceleration and emission models 

Non uniform velocity fields in collisionless plasmas allow 
a variety of particle acceleration processes where the density 
in high energy particles can grow until it affects the flow it- 
self and the acceleration mechanism. Indeed, most of the 
VHE sources detected so far harbor powerful flows, especially 
AGN which appear as perfect laboratories to study relativis- 
tic plasmas under extreme conditions inaccessible to experi- 
ments. Fermi acceleration processes of the 1st and 2nd or- 
der in shocks and turbulence, respectively, have been exten- 
sively explored in the literature and are clearly expected in 
AGN lEHJ |99l [lOOl HOH HSU UnS] UMl USSll, as well as mag- 
netic reconnection phenomena 11061 11071 . and possibly di- 
rect electric fields in gaps or centrifugal forces around rotat- 
ing SMBH (89J [JU81 [J091 [JTUl [TTTJ [TT2) . There is a growing 
interest especially in relativistic reconnection in recent years 
since intermittent reconnection events in sheets or small mag- 
netic islands in jets or in the SMBH vicinity could be recog- 
nized as efficient particles accelerators and nicely match AGN 
non-thermal emission properties, especially fast VHE variabil- 
ity m^rTT^nT5irTT6inT^riT^nT^n^oin^Tin^arT2^n^4i . 

Together with cosmic ray and neutrino experiments, CTA is 
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among the most promising projects for detailed studies of these 
extreme cosmic accelerators. 

VHE observations provide an uncluttered view of the physics 
by pointing out only the most extreme phenomena at the highest 
energies and by disentangling non-thermal effects from ther- 
mal ones radiating at longer wavelengths. They should help 
identifying dominant processes in the global energetics of the 
source. Powerful phenomena inducing VHE emission are po- 
tentially the sources of several other secondary events observed 
at longer wavelengths, the explanation of which possibly lies in 
the primary VHE data. CTA will therefore offer an invaluable 
tool for exploring the violent mechanisms at work in SMBH 
environment, especially on the jet physics, including forma- 
tion, collimation, and propagation as well as the development 
of shocks and turbulence along the jet. Depending on the pre- 
cise origin of the VHE radiation, which is still to be settled, 
CTA data could also constrain properties of the SMBH and the 
accretion regime. 

As presented in part 2, several scenarios have been proposed 
to explain the TeV emission of blazars but none of them is yet 
fully self-consistent. In the absence of a definitely convincing 
global picture, a first goal for CTA will be to constrain model- 
dependent parameters in a given scenario. For instance, basic 
stationary one-zone SSC scenarios have typically eight main 
free parameters, three macroscopic ones for the size, the mag- 
netic field, and the Doppler factor of the emitting zone, and five 
others to describe the microphysics of the emitting particle dis- 
tribution (density factor, slopes of a broken power-law, maximal 
individual Lorentz factor, and Lorentz factor at the break). High 
quality VHE spectra and detailed simultaneous SED are needed 
during quiescent as well as active states to get a similar num- 
ber of observational constraints (as frequencies and fluxes at the 
synchrotron and the inverse-Compton peaks, break frequencies, 
and spectral indices). This will be achievable thanks to the large 
spectral range, high sensitivity and high spectral resolution of 
CTA, together with the coordination of multi wavelength cam- 
paigns. In such a way, the physics of basic radiation models will 
be well constrained by CTA, characterizing the properties of the 
particle distribution from its emission at VHE, which should 
help to identify the acceleration mechanisms at work, and de- 
cide which part of the models should be corroborated or ruled 
out. 

The second most difficult goal will be to distinguish between 
the different remaining options and to firmly identify the dom- 
inant acceleration and radiation mechanisms. Detection of any 
specific spectral features, break, cut-off, absorption or addi- 
tional components which could be radiative signatures of spe- 
cific acceleration processes, would be mandatory in this regard. 
Observations above several TeV will be an important test for 
the most popular SSC scenarios which may have difficulties to 
generate hard enough spectra because of the limitation due to 
the Klein-Nishina regime. This could rule out simple one-zone 
models by requiring alternate photon sources from, for exam- 
ple, a structured jet. Such intense research will benefit from 
improved performances of current particle-in-cell (PIC) codes. 
Probing AGN variability at all time scales, down to the shortest 
ones at the few-seconds scale with CTA, will significantly con- 



strain acceleration and cooling times, instability growth rates, 
shocks and turbulence time evolution. Deep analysis of the light 
curves should explore all timescales from several years down to 
several seconds, look for periodic and quasi-periodic behaviour, 
and characterize the additive or multiplicative processes trying 
to distinguish macrophysics from microphysics effects at the 
origin of the variability. Periodicities of about 10 to 50 min- 
utes, related to the last stable orbit around SMBH of 10 7 so- 
lar masses, should be reachable with a good temporal coverage 
with CTA even during low activity states (see section 5). The 
role of CTA as a time explorer will be decisive for constraining 
both radiative phenomena, and global geometry and dynamics 
of the AGN engine. 

4.1. Testing and constraining leptonic scenarios with CTA 

In this section, we illustrate some aspects of leptonic scenar- 
ios that the unique capabilities of CTA should be able to char- 
acterize, confirm or invalidate. Leptonic models assume the 
emission of relativistic electrons with energy E = ym e c 2 and 
individual Lorentz factors up to about 10 6 . First order Fermi 
acceleration processes in shocks can for instance convert about 
10% of the bulk jet kinetic energy into random energy of fast 
particles and provide a power-law particle energy distribution 
N(y) oc y~ n required to explain the observed power-law spectra 
of blazars, with a synchrotron flux F sya oc v~ a , with spectral 
index a - (n - I) 12. These electrons can up-scatter their own 
synchrotron photons (SSC) or other externally produced low- 
frequency photons (EIC), increasing their energy by several or- 
ders of magnitude and inducing strong VHE radiation. 

4.7.7. Peak frequency and emission level correlation in SSC 
models 

A noticeable characteristic of SSC scenarios resides in the 
correlation which should be found in the evolution of the high 
energy IC peak. Indeed, such types of correlation have been 
already found at lower energies for the synchrotron bump, as 
shown for instance by a detailed analysis of the UV and X- 
ray spectra of the well-known BL Lac Mrk 421 012511 which 
display evidence for a clear correlation between the position 
of the synchrotron peak (£ pe ak) and the emission level at the 
peak £peak « v^™2 * . However, such effects remain out of 
reach for many sources at VHE with current IACT because their 
study requires precise monitoring of the spectrum over a wide 
range of energies. This is especially difficult in the case of high 
energy peaked sources where the IC peak is located close to 
the high energy threshold of orbital observatories (e.g. CGRO, 
Fermi) and low energy threshold of ground based Cherenkov 
instruments (e.g. H.E.S.S., MAGIC, Veritas). 

CTA with its improved low energy threshold and high sen- 
sitivity will offer the unique opportunity for detailed investiga- 
tions of the IC peak evolution of bright VHE blazars such as 
Mrk 421, Mrk 501 and PKS 2155-304. Detecting the peak of 
the VHE bumps and monitoring its temporal evolution on short 
timescales for a small sample of blazars will be a real bench- 
mark to constrain emission models, which can be achieved by 
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CTA. Improved CTA performances, or the hypothetical obser- 
vation of an exceptionally bright event, could extend such anal- 
ysis to the family of radiogalaxies, which would be a significant 
advance. 

Indeed, high quality data on the evolution of the synchrotron 
and IC peaks can drastically constrain SSC modelling. As 



injection & cooling 



shown in Fig 12 the short timescale evolution of the IC peak 



during flares is directly related to the physical origin of the vari- 
able event and can characterize it. In an "injection & cooling" 
scenario, the evolution of the particle energy spectrum (N(y, t)) 
can be described by the kinetic equation 



dmy, t) 

dt 



By 



[{%yn + Yic}N(y,t)] = Q(y,t) 



(6) 



where y syn and jiq are the synchrotron and IC cooling rates re- 
spectively. Here non-radiative losses and escape are neglected. 
The particle acceleration is a free term, just described by a given 
injection rate Q(y, t). In such a case, the density of high en- 
ergy particles significantly increases inside the source during 
the injection phase, with a quick increase of the observed emis- 
sion. When the acceleration becomes inefficient and the injec- 
tion stops, the radiative cooling takes control over the particle 
energy evolution and reduces the number of high energy par- 
ticles. The two processes, "injection & cooling", can explain 
observed flares giving characteristic "1" like shape for the peak 



evolution, as shown in Fig. 12 (top panel). This model predicts 



that the correlation should be different for the rise and decay of 
a flare. This can characterize some AGN events, but not all of 
them, remembering the case of Mrk 421 where the correlation 
measured for the synchrotron bump M125I suggests that they are 
similar. 

In an "acceleration & cooling" scenario, one can assume that 
the particles are systematically accelerated inside the source, 
following a kinetic equation of the form 

dN{yJ) + l[[A(y,t)-(y^ + y lc )}my,t) 



dt 



dy 



+D(y,t) 



dN(y, t) 



dy 



= 0, (7) 



where two additional terms describe acceleration (A(y, t)) and 
stochastic diffusion (D{y,t)) in the energy domain [126]. The 
maximum energy particles can reach depends on the efficiency 
of acceleration versus radiative cooling and on the duration of 
the acceleration process. When acceleration becomes ineffi- 
cient, cooling starts reducing the particle's energy. The evo- 
lution of the particle energy spectrum then draws a "wing" like 



shape for the IC peak path, as shown in Fig. 12 (middle panel) 



In this case as well the correlation should be different for the 
rise and the decay, except for a fine tuning of the parameters. 

In the alternative case where the variable event is simply due 
to a change of the source Doppler factor, and not to the evo- 
lution of the particle energy distribution, the same peak cor- 
relation occurs for the rise and decay of a flare, as shown in 
Fig. 12 (bottom panel). Such scenario can be relevant when the 



emitting zone travels on slightly curved trajectory and its ap- 
parent Doppler factor 6 significantly varies for different values 






^ 

*- 



i 



10.0 
10.5 
11.0 
11.5 
12.0 

12.;, 



-13, 







c 

- 
— 


~]5^v\ - 


l^—^^^ 




~T~\\ \\ " 


— CTA-B, f,=15 min 















a - 






:4.0 24.5 25.0 25.5 26.0 26.5 27.0 27.5 
lo Kl3 IVi [Hz] 

acceleration & cooling 



-Hi 









— 


CTA-B, 6—15 min 


r***. 


||p\ 




/\ \ \l- 




i \ \ i 1 ii\ 1 



22 



2'1 



24 2. r > 26 

logioM (Hz] 
beaming Qffgcl 



27 



2ii 




09 



24 ?5 26 

logiuW [Hz] 

Figure 12: Examples of the IC bump and peak evolution for various time- 
dependent SSC scenarios. Spectra during rise (resp., decay) time are shown 
at different moments in red (resp., blue) thin lines. The IC peaks of the spectra 
are connected with red and blue bold lines which show a shape characteristic of 
the different scenarios, namely a '"1" like shape characteristic of an "injection 
& cooling" scenario (top panel), a "wing" like shape line drawn by the IC peaks 
for an "acceleration & cooling" scenario (middle panel) and a "single line" for 
a scenario with a simple change of the Doppler factor (bottom panel). The three 
examples given here represent typical bright VHE flares expected from nearby 
TeV BL Lac sources. The middle panel specifically reproduces the rise and 
decay of a Mrk 501 flare 1 127] while the bottom one describes an active event 
of Mrk 421 over a few hours 1 18 1. Bold black lines show the CTA sensitivity 
curve for 15 minutes of integration time with the B array. Only the high sen- 
sitivity of CTA provides temporal resolution good enough to follow such fast 
sub-hour evolution. 



of the viewing angle 6. Since the observed flux is proportional 
to a high power of the Doppler factor a small change of 9 (e.g., 
from 5 to 3 degrees) induces a significant increase of the ob- 
served emission 111281 . Current available data do not constrain 
well the evolution of the VHE bump since the IC peak cannot be 
monitored yet, being below the threshold of current IACT and 
difficult to measure on short time scale by Fermi. As clearly 
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Figure 13: Simulations of spectra obtained in 15 minutes by CTA (B array) 
during the rise of the flare for the case of the simple change of Doppler factor 
illustrated in Fig |12| (bottom panel). This illustrates the decisive importance of 
having good sensitivity and spectral resolution over a large range in frequencies, 
especially at low energies, in order to allow a good determination of the peak 
motion during the flare. 



shown in Fig. 13 the performances aimed with CTA should 
allow for the first time to follow these behaviors. 

4.7.2. TeV and X-ray emission correlation 

One strong piece of evidence in favour of SSC models is the 
correlation often found between the evolution of X-ray and the 
VHE flares. This is not necessarily the case for EIC scenar- 
ios where there may be a competition between the two types 
of radiative losses, synchrotron and external Inverse-Compton, 
which can even result in an anticorrelated behavior. However 
some hadronic models, such as the synchro-proton scenario 
where protons and electrons are simultaneously accelerated, 
can also explain such a correlation, though this deserves a more 
detailed investigation. 

This correlation is not universal, as definite orphan VHE 
flares have been observed, with no X-ray counterpart (and con- 
versely). Light curves in X-rays and gamma-rays can some- 
times look very different, as shown for instance by Fermi for 
3C279 01291 . In this source, comptonization of external radi- 
ation from the dusty torus and the broad line region should be 
taken into account for the generation of gamma-rays, and X- 
rays seem to have a different origin than in HBL, the steady 
X-ray emission being possibly an inverse Compton component 
corresponding to a synchrotron bump peaking in the mm/sub- 
mm range. Intensive and long term coordinated monitoring in 
VHE and lower frequencies should throw more light on these 
questions, studying the variation of the SED as a function of 
source activity [130]. Pursuing X-ray campaigns will be essen- 
tial in this regard, and will contribute to the characterization of 
the emitting particle distribution. 

4.1.3. Hard intrinsic spectra 

Despite the present uncertainties on the EBL values, some 
AGN observations suggest that intrinsic TeV spectra could be 



unexpectedly hard, with intrinsic energy spectral index a em 
equal to 0.5 or even smaller for mid or high levels of the EBL 

mm 

For instance, the observations of 1ES 1101-232 made by 
H.E.S.S. show a TeV spectrum that can be approximated by 
a power-law function with spectral index a b s = 1.88 ± 0.17. 
The emission of this distant (z - 0.186) source is absorbed by 
the EBL and current uncertainties on the EBL flux level and 
spectrum 11321 11331 1134 1 introduce some difficulties in defin- 
ing how hard the intrinsic spectrum is. However, in the case of 
1 101-232, it tends to be very hard, with intrinsic spectral index 
a ~ 0.5 when corrected for low EBL levels, and higher than 
0.5 if the true EBL was significantly higher than its strict lower 
limit derived from galaxy count [ 131111351 . Such hard intrinsic 
spectra could be also expected in a number of extreme blazars 
like 1ES 0229+200 and 1ES 0414+009 fJ9l[T36l . 

At a first glance, it seems difficult to generate such hard spec- 
tra within simple SSC models and standard particle acceleration 
models. Then, the IC scattering that gives the TeV emission 
usually occurs in the Klein-Nishina regime and the efficiency 
of the scattering decreases very quickly with increasing elec- 
tron and photon energies, resulting in a soft intrinsic spectrum 
witha em < 1)- 

A simple solution to this problem is to assume a sharp 
low energy cut-off in the particle energy spectrum |137]. If 
the minimum energy of the particles is high enough (E m \ n - 
7min'MeC 2 ,7 m i n > 10 5 ) then the intrinsic spectrum can be hard 
with spectral index o- em > -1/3. This limiting value of the 
spectral index comes from the fact that the particles are scatter- 
ing 'tail' photons of the synchrotron emission. In such approach 
the value of a em depends on the y m j n value and can extend from 
-1/3 up to 1. Moreover, the scenario predicts an abrupt break in 
the spectrum at energies of several TeVs. Therefore, it should 
be easy to test it, using CTA to obtain large band spectra from 
the GeV range up to a few tens of TeVs. Various ways to gener- 
ate such specific particle distributions and hard VHE spectra are 
presented in [39|. Such investigation aims to find a simple and 
robust description for the VHE emission of the TeV blazars, to 
be able to predict their intrinsic spectra, to relate them to the 
AGN characteristics, and to better constrain intervening media. 



4.1.4. Time delay in TeV light curves 

Detailed observations of the rapid activity of Mrk 501 con- 
ducted by the MAGIC telescope on July 9th, 2005 Ifl38l show 
delays between the light curves obtained in different energy 
ranges. In particular a several minutes delay was seen between 
the light curve obtained in the range 150-250 GeV and the one 
in the range 1.2-10 TeV. A natural explanation for such delay is 
an on-going particle acceleration process just caught at the time 
of the observation, as proposed by 1139111401 . Quantum nature 
of the gravity is an alternative possibility [141]. Such type of 
delays deserves precise investigation II142II1431I1441 and more 
advanced data, especially the spectral shape evolution is crucial 
to confirm or reject the acceleration scenario. This will be pos- 
sible only with CTA, with sensitivity good enough to provide 
high quality spectra of such rapid flaring events. 
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4.2. Hadronic and lepto-hadronic models 

Observational evidence and theoretical considerations on the 
energetics of AGN 11451 11461 [T47l suggest that AGN jets can 
contain a non-negligible fraction of hadrons which can consti- 
tute an important fraction of the particle content. The character- 
istics of these hadrons, e.g., their energy spectrum and number 
density, are however not known and it is not yet clear if they 
contribute significantly to the radiative emission. 

Hadronic emission models explore scenarios where radiative 
processes involving relativistic hadrons contribute to the high 
energy emission from AGN. These scenarios provide a direct 
connection between the detectable gamma-ray emission and the 
(as of yet undetected) emission of neutrinos and ultra-high en- 
ergy cosmic rays (UHECRs) from AGN. Despite difficulties of 
establishing a firm link with the arrival direction of UHECRs, 
AGN remain one of the few candidate sources for these parti- 
cles, and the nearby radio-galaxy Centaurus A is still consid- 
ered a potential source for several UHECR detected with the 
PAO [ 148 ]. Current and future neutrino telescopes, such as Ice- 
Cube and KM3NeT, might detect very energetic neutrinos from 
certain classes of AGN, which would provide firm evidence that 
hadronic processes take place in these objects. The higher tem- 
poral and energy resolution of CTA and its improved sensitivity 
will help investigate such scenarios. 

4.2.1. Hadronic emission scenarios 

Different hadronic models have been developed to describe 
the SEDs of different types of blazars and radiogalaxies. They 
usually require magnetic fields of the order of 10 to 100 G, 
much higher than in the leptonic description. The Synchrotron 
Proton Blazar (SPB) model Q49] [BO] [BJ] [B2 ascribes the 
VHE emission in HBLs mainly to synchrotron emission from 
ultra-relativistic protons, pions and muons, while the X-ray 
emission is dominated by electron synchrotron photons. The 
"Proton Induced Cascade" models M153II1541 . in which photo- 
meson (mainly pion) production by ultra-relativistic protons 
with the internal radiation field in the source dominates the 
VHE peak, are better adapted to describe the SEDs of LBLs 
than that of HBLs. In FSRQs, photo-meson production on ex- 
ternal photons from the disk, the BLR and the dust torus adds 
to the other hadronic emission processes [155]. 

A general problem for hadronic models is that short-term 
variability appears more difficult to account for in hadronic than 
in leptonic scenarios given the anticipated longer time scale for 
the interaction processes. However, if the variability is orig- 
inated by changes in the viewing angle of the jet, the effect 
due to the variation of the bulk Doppler factor is the same for 
any process taking place in the jet's frame Ml 561 11571 . They 
need also some elaborated solutions to explain the TeV and X- 
ray emission correlation, which is not expected in the simple 
case of co-acceleration of electrons and protons in the same 
region. Indeed, alternative scenarios have been proposed that 
might be able to account for rapid variability and correlated 
flares (BJ [BJ [BJ [BJ] [BJ [BJ] [BJ [B2. One recent 
example [164] in which AGN jets interact with gas envelopes 
from red giant stars and emit proton synchrotron radiation can 



even explain the very rapid variability seen during the big flare 
ofPKS 2155-304 in 2006. 

4.2.2. Testing hadronic scenarios with CTA 

Although it is not obvious to find a clear detectable signature 
that would allow the unambiguous distinction between leptonic 
and hadronic models, the highly improved information from 
CTA on the SED and its evolution, together with multiwave- 
length and multi-messenger data from other instruments, will 
help constrain model parameters and reject certain scenarios. 
An interesting target for such studies would be sources with 
significant internal or external photon fields, such as LBLs or 
FSRQs, in which photo-meson production and cascading play 
an important role in hadronic scenarios. Figure[T4](upper panel) 
shows one scenario, where cascades lead to a hardening of the 
SED at the highest energies, which will be clearly detectable 
with the better energy coverage and sensitivity of CTA. The 
search for spectral variability could provide another piece of 
information for the existence of cascading. 

Photo-meson production in internal or external photon fields 
would lead to a characteristic cutoff energy of the proton spec- 
trum, depending on the photon field density. This cutoff could 
be inferred from observations of the proton-synchrotron emis- 
sion. Given the large number of sources accessible to CTA, it 
would be possible to study a relation between the cutoff energy 
and the blazar class, thus probing the blazar sequence within 
the hadronic framework. A refined understanding of the EBL, 
which could also induce features at the high energy end of the 
VHE spectrum, will be required and should be developed in 
parallel with the studies proposed here. 

Apart from the two extreme cases, where the emission in the 
high energy bump stems either from the leptonic or from the 
hadronic population, one can also investigate mixed scenarios, 
in which both leptonic and hadronic processes contribute to the 
high energy bump. Such "lepto-hadronic" models can exhibit 
detectable spectral features in the high energy and VHE range 
as a result of the combination of different components. An il- 



lustrative example is shown in Figure 14 (bottom panel). These 
different components would also be expected to exhibit differ- 
ent temporal behaviour that one could hope to resolve with the 
high timing capabilities of CTA. 

For all these studies, simultaneous multiwavelength coverage 
is needed to assess the electron population in the source from 
observations of the synchrotron bump. With this information, 
the existence of an additional hadronic component can be eval- 
uated from an investigation of the VHE emission with CTA. 

4.3. Studying spectral features 

In a few cases, VHE spectra are already known to deviate 
from a pure power-law. This is better seen for bright sources 
during active states. Spectral break and curvature have been 
detected in the SEDs of Mkn 421, Mkn 501 and PKS 2155-304 
II167II1681I169I . The origin of such spectral features is not yet 
well identified and will require high quality spectra obtained 
with CTA. They could be the signature of different phenomena 
such as a cut-off due to the maximal Lorentz factor y max of the 
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Figure 14: Reproducing SED with hadronic and lepto-hadronic models. Up- 
per panel: The SED of the FSRQ 3C 279 compared to a hadronic model 
with several individual radiation components, namely proton synchrotron and 
cascade (dashed), muon-synchrotron and cascade (triple-dot-dashed), cascade 
from neutral pions (dotted), and cascade from charged pions (dot-dashed), 
from [152) (see also (165)). Bottom panel: The SED of the HBL PKS 2155-304 
as seen during the multiwavelength campaign of 2008 [30], compared to a lep- 
tonic model (red line and dotted green lines for the synchrotron and SSC com- 
ponents). An additional component from proton-induced cascades is shown 
in magenta (dashed dotted lines) 1 166|. High sensitivity and spectral resolution 
are needed to look for signatures of the various hadronic components with CTA. 



population of emitting particles, to the Klein-Nishina regime, or 
to some absorption effects. The detection of additional bumps 
in the SED could provide evidence for the presence of some 



hadronic component (see Fig 14 1. 

Indeed, recent observations by the Fermi Gamma-ray Space 
Telescope (Fermi) Large Area Telescope (LAT) revealed that 
the spectra of bright blazars cannot be described by a sim- 
ple power law model M170I 11711 . A much better description 
is obtained with a broken power law, with break energies of 
a few GeV. In the brightest GeV blazars, the photon statistics 
allows spectral variations at different flux levels to be stud- 
ied. In spite of obvious spectral changes, the break energy in 
some sources like 3C454.3 turned out to be surprisingly stable 
II172II173TI174I . These breaks can be produced by absorption of 
the GeV photons by the photon-photon (yy) pair production on 



BLR Q75J (Fig. [T5| 
energies 261 GeV/(£ , LyC (eV)[l +z]), i.e. at about 5/(1 + z) and 
19/(1 + z) GeV. It has been argued that a relatively high opacity 
in the He n LyC observed in a number of bright blazars implies 
the location of the y-ray emitting region within the highly ion- 
ized inner part of the BLR B175I . The absorption internal to the 
BLR thus complicates the determination of the intrinsic spectral 
shape of the bright blazars. The broad-band coverage provided 
by Fermi and CTA will be crucial for such a study. 

Photoionization models of the BLR predict also other strong 
lines, with Balmer lines (Ha, H/3) and Paschen a being the 
strongest. These lines can produce absorption breaks at 100- 
150/(1 + z) and 400/(1 + z) GeV, respectively B176I . i.e. within 
the CTA range. Thus with a spectral resolution of about 10%, 
the breaks can easily be detected. The detailed fitting of the 
broad-band spectra including internal absorption within the 
BLR will give information on the column density of various 
soft photon sources that produce absorption in a given gamma- 
ray band. Thus, CTA will be able to map the BLR structure 
and to locate the gamma-ray emitting region with a higher pre- 
cision than was possible up to now. Highly active states should 
provide the opportunity to obtain significant detections of such 
spectral features. Such type of studies will clearly be very sen- 
sitive to the low energy threshold of CTA. 

5. The origin of variability 

5.1. Timing capabilities of CTA 

Blazars are usually detected in the gamma-ray band when 
in a high state, though the sensitivity of current IACT has re- 
cently allowed the community to start probing the quiescent 
state of some of the strongest VHE sources 111691 . The high 
states of these objects are dominated by strong and burst-like 
flaring episodes that are characterized by short, aperiodic vari- 
ability events possibly dictated by the fast cooling times of the 
~ 10 TeV electrons thought responsible for the VHE emission. 
The detection and detailed studies of these short transient events 
is important to put constraints on the physics of the emission 
process and the astrophysics of AGN jets and is one of the key 
science cases for CTA. 

Its huge effective area (10 5 to 10 7 m 2 , much larger than the 
m 2 effective areas reachable by space experiments) makes CTA 
an ideal instrument for the study of timing properties of blazars. 
Detections of fast variability of the VHE emission from some 
AGN indicate that their emission region is extremely compact. 
Such a compact emission region would most naturally occur 
close to the SMBH, in the blazar central engine. Then timing 
properties of the AGN flares encode valuable information about 
the properties of the black holes and their immediate surround- 
ings. However fast variability could also indicate that there are 
very compact substructures within the jet. Current generation 
instruments are able to detect variability on the scales of sev- 
eral minutes in the case of the brightest blazar flares. Detect- 
ing characteristic temporal features with CTA, such as minimal 
variability times, or periodic or quasiperiodic events, could pro- 
vide definite answers to the question of the size and possible 
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Figure 15: Top: Spectral energy distribution of a few blazars as observed with 
Fermi/LAT during the first 6 months of its operation. The best-fit broken power 
law and a power law with the double-absorber models are shown by the dashed 
and solid lines, respectively. Bottom: Spectral energy distribution of 3C 454.3 
at low and high fluxes as well as averaged over the whole observation period 
of 2.5 years. The solid curves show the best-fit lognormal model with absorp- 
tion by He ii and H Ly recombination continua. The dotted curves show the 
unabsorbed lognormal distributions. 



location of the VHE emission region (active zone in extended 
or VLBI jets, base of the jet, or central engine magnetosphere) 
and stringent information on Doppler factors, physical emission 
mechanisms, particle acceleration, cooling and escape times, as 
well as on the intrinsic time scale of the central SMBH. 

Variability time scales of AGN time series provide direct 
constraints on the y-ray emitting zone. Short rise/decay times 
of an observed burst set an upper limit on the minimum vari- 
ability time of the source, constraining by causality its size 
and Doppler factor. Complementary to variability studies in 
the temporal domain, one can derive the Fourier properties of a 
time series e.g., from Power Spectral Density (PSD) in a given 
frequency range, set by the observation live time and the sam- 
pling frequency of the lightcurve. Such analyses character- 
ize the stochastic mechanism generating the flux modulation. 
These studies, inspired by the X-ray domain [177|, are now 
feasible in the TeV range for a handful of dramatic outbursts, 
such as the exceptional flares of the blazar PKS 2155-304 seen 
by H.E.S.S. in July 2006 [2| which produced the fastest and 
highest fluence flares ever observed in the VHE range from an 
AGN, and can thus serve as a proxy to study what are the lim- 
iting variability scales that CTA is likely to be able to probe 
from these kinds of objects. Nonetheless, the lack of sensitiv- 
ity and limited field of view of the current generation of IACT 
are responsible for their scarcity. The next generation, repre- 
sented by CTA, will allow for faster sampling of the lightcurves 
and for the generalization of studies in both Fourier and tempo- 
ral spaces. The expansion of this field will be enabled by the 
increased sensitivity and the lower energy threshold that CTA 
will afford (eg. ~ 200 GeV for H.E.S.S. compared to 50 GeV 
for CTA). Hereafter we illustrate the consequences of these in- 
creased capabilities by simulating PKS 2155-304 flares as pos- 
sibly monitored with CTA. 

The lowering of the energy threshold will indeed allow one 
to monitor the source between 50 GeV and 200 GeV, whose 
flux is modeled here using the spectral energy density derived 
during the 2008 multiwavelength campaign on PKS 2155-304 
1130111781 . The increased energy coverage, combined with the 
large simulated effective areaj lfT31 11791 , will allow CTA to 
significantly detect the source within a shorter duration than 
any current IACT, hence extending the high frequency part 
of the PSD. The low frequency part of the PKS 2155-304 
spectrum during the high state is well described by a power 
law of index 2 1 169 180|. Additional variance in the lightcurve 
beyond the reach of H.E.S.S. is derived from the Timmer and 
Konig method 11811 . assuming the PSD extends to higher 
frequencies with the same index. Finally, the binning of the 
simulated CTA lightcurve is performed assuming an average 
significance that equals the one of the monitored H.E.S.S. 
lightcurve. 



One specific CTA lightcurve realization is shown in Fig. 16 
Provided the red noise behavior extends to higher frequencies, 
additional substructures which could not have been resolved 



" The final results are barely dependent on the configuration used, yielding 
an uncertainty on the final temporal binning of around 5%. 



17 



by H.E.S.S. are revealed by the peak finding and fitting 
procedure. During such an event, the shortest significant 
rise time accessible to CTA would be r r = 25 + 4 s (sec- 
ond peak on Fig. 16 1, which is approximately seven times 
smaller than r r h.e.s.s. = 173 ± 28 s. Such an upper limit 
on the smallest variability time scale constrains the ratio of 
the size of the emission region over the Doppler factor to 
R6~ l < cr r /(l + Z) ~ 6.7 x 10 n cm, implying e.g. an emission 
region smaller than 1 AU for a Doppler factor of 20, an unusual 
value in the current acceleration models used for active galactic 
nuclei. The detection of such short time scale events would be 
a real challenge for jet formation models. 




Figure 16: Simulated integral flux of PKS 2155-304 above 50 GeV as CTA 
would monitor it. This simulation relies on an extension of the red noise behav- 
ior to high frequencies, generating the short time scale structures (second and 
fourth peaks). The data are binned in 7.5 seconds intervals. 



This first example illustrates the great CTA timing capabil- 
ities which could detect events almost ten times shorter than 
the shortest ones detected by current IACT This temporal 
resolution, below the minute time scale under exceptional 
circumstances, could impose severe constraints on acceleration 
mechanisms and raise the question of the maximal Fourier 
frequency accessible to very high energy sources. However, 
such events are relatively rare and will require long term 
monitoring of large samples of AGN and the triggering of 
observations under alarm to be efficiently detected. Below we 
present a more general statistical study which aims at assessing 
the typical gain in terms of temporal resolution of flaring events 
that CTA could monitor. 



mization of the Bayesian priors of the analysis is described in 
the Appendix B. 

Two groups of 1,000 light-curves were derived for H.E.S.S. 
and for the CTA extrapolation. Observe that the power law de- 
scribing the PSD being essentially the same for both groups, 
the only difference between the two data sets is the sampling 
provided by the instrument. The distribution of flare doubling 
times tj for both samples is shown in Fig[l7] Many of the 
light-curve features not resolved in the H.E.S.S. sample can be 
separated in different events with the CTA resolution, shifting 
the peak of the distribution towards shorter duration events and 
thus revealing additional information about the timing structure 
of the source which was previously unnoticed. In particular, the 
most probable flare r^ derived for H.E.S.S. is of the order of 5 
min, a typical value for the H.E.S.S. light-curve of PKS 2155- 
304 in 2006. With CTA, this value is expected to shift to 100 s, 
a factor of 3 shorter than H.E.S.S. 
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Figure 1 7 : Kernel density estimates for the distributions of the sizes of variabil- 
ity features detected from the H.E.S.S. and CTA simulated light-curve samples. 
The most probable typical flare size as seen with CTA, of the order of 100 sec- 
onds, is approximately l<x from the typical values of 5 minutes observed today. 



To perform this statistical study, we chose to adopt un- 
binned, non-parametric methods, in order to take into account 
all the available timing information present in the light-curves 
and allow to probe for short variability events in a model- 
independent way. We generated a family of light-curves [ 181 1, 
with an empirical PSD index /J ~ 2, which are independent 
realizations of the process describing the H.E.S.S. observa- 
tions of 2006. A second group of light-curves was also gen- 
erated, extending the PSD to the higher frequencies which 
would be accessible by CTA, as described in our first example. 
These unbinned sequences of time-tagged events were further 
analyzed for variability using the Bayesian Blocks algorithm 
|[182l [1831 [1841 [1851 . which aim at finding the best partition 
of a sequence of events in a non-parametric fashion. The opti- 



The probability of occurrence and detection of a given ultra- 
fast flare as presented in the beginning of this section can be 
assessed by comparing the cumulative distribution of the flare 
durations as derived from our MC simulations (Fig[T8J). If short 
events are present in PKS 2155-like light-curves, these will be 
better sampled by CTA, up to the high resolution limit of the in- 
strument. In the PKS 2155-304 light-curve as seen by H.E.S.S. 
in 2006 [2 1, the shortest variability event was of 173 s, which 
from simulations has an estimated probability of occurrence 
and detection by an instrument like H.E.S.S. of ~ 13% (a rare 
event). With the improved capabilities of CTA, such an event 
will be rather usual, with a probability of ~ 40%. The short- 



est burst shown in Fig 16 as simulated for CTA, could not be 
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Figure 18: Cumulative function of the same distributions shown in Fig |17| (see 
text). 



detected by H.E.S.S. due to the limitations in sensitivity of the 
instrument. CTA, though, will be capable of seeing such a 25 s 
event if it is expected to be a rather rare occurrence, with prob- 
ability below 10%. 

To conclude, the advanced timing performance of CTA will 
allow to sample AGN light curves up to timescales significantly 
shorter than the ones achievable by the current generation of in- 
struments, with the potential to probe variability well below the 
minute-scale for exceptional outbursts. The statistical study of 
typical simulated AGN light-curves shows that CTA will de- 
tect routinely many more variable events than current IACT, as 
a direct consequence of the better sensitivity, which will allow 
to resolve previously unnoticed short-time variability features. 
Typical flares duration are expected to be about 3 times shorter 
than those of the current generation, reflecting better the in- 
trinsic structure of the sources. One can anticipate that such 
improvements in the timing analysis capabilities with the next 
generation of IACT observatories will be of great relevance for 
the study of AGN jet physics, bringing us closer to the fun- 
damental limit of these variability timescales and to an under- 
standing of the physical processes ultimately governing them. 

5.2. Multiwavelength and multi-messenger campaigns: moni- 
toring and targets of opportunity (ToO) 

Given the complex variability properties of AGN, the flexi- 
bility of the CTA Observatory will be extremely useful. Long 
term correlated multiwavelength monitoring programs of a 
well-defined sample of AGN will be a necessity for several 
reasons, first to characterize quiescent and low-activity states 
and to estimate their VHE duty cycle, and also to send alarms 
for ToO during specific events (bright flares, rare events seen 



at lower energies or with other messengers, and also particu- 
larly low-activity states for instance in the optical range, which 
should be a good time to determine redshifts of remote BL Lac 
objects when spectral features from the host galaxy can be bet- 
ter detected above a low non-thermal optical continuum). The 
possibility to observe with CTA sub-arrays will be mandatory 
for such observing programs 1 186]. 

Characterizing the VHE variability properties themselves 
should help identifying the mechanisms at the origin of the tem- 
poral evolution, and to distinguish between features induced 
by specific bursts from those essentially due to a noise pro- 
cess. Multiwavelength monitoring should also contribute to 
identify the various emission zones, to see whether they can 
be co-spatial or not, to locate them relatively to the others. This 
has been recently illustrated 11871 11881 . combining the GeV 
gamma-ray lightcurves from Fermi satellite with X-ray, optical 
and VLBA monitoring data. With current generation instru- 
ments the VHE gamma-ray long term monitoring is possible 
only for handful of brightest AGN. With the CTA instead we 
should be able to clarify whether the VHE variability is domi- 
nated by the jet physics, with development of shocks and turbu- 
lence in the jet itself, or if the origin of the variability primarily 
lies in the accretion flow and the immediate vicinity of the cen- 
tral black hole, and to disentangle variability due to the global 
dynamics of the engine from variability induced by radiative 
processes (see section 4) studying the variation of the SED as 
a function of source activity. As an example, the binary black 
hole system OJ 287 could be a unique laboratory in this regard, 
thanks to the regular modulation of its parameters due to the 
orbital motion, which allows to observe the system, complex 
but rather well constrained, in different conditions over years 
[189|. Long timescale observations in X-rays and multiwave- 
length succeeded to answer such questions in the case of Seyfert 
galaxies |190|. Establishing for instance the relationship be- 
tween the RMS variability and fluxes, they ruled out simple 
noise models and showed the presence of non-linear variability, 
possibly driven by fluctuation in the accretion flow inducing 
a multiplicative process 11911 . Multiplicative processes have 
recently be suspected in VHE blazars, as shown by temporal 
analyzes of the big flare of PKS 2155-304 |fT69l . Variability 
studies of TeV AGN are still in their infancy, with variability 
characteristics still poorly determined and difficult to interpret 
[ 192|. As shown above, CTA will open an avenue to deep time 
analysis at VHE. 

Generally speaking, AGN being very broad-band emitters 
over the whole electromagnetic spectrum, very promising syn- 
ergies can be anticipated with most space and ground-based 
large infrastructures of the domain such as LOFAR, SKA, 
GAIA, JWST, and ELT Based on the results obtained so far 
by the MOJAVE project and similar VLBI searches |[T93l[T94l 
195|, one might expect that regular VLBI mapping at radio fre- 
quencies, coordinated with VHE campaigns, will significantly 
constrain the global picture of the AGN central core at the sub- 
pc scale, linking the mechanism of jet formation, the ejection 
of radio knots, the growth of shocks and turbulence together 
with the localization, geometry and dynamics of the VHE emit- 
ting zones. It should also solve the persistent paradox of non- 
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detection of superluminal velocities in the VLBI maps of VHE 
blazars while extremely large bulk Lorentz factors are needed 
for explaining the fluxes and variability [196]]. In the optical 
range, monitoring of the photometric and polarized continuum 
fluxes, and weak spectral lines from the broad line regions (if 
any), will allow to analyze the unclear link between VHE and 
optical active states, to test the relative importance of EC contri- 
bution in the blazar sequence and especially FSRQ, and to make 
the connection between accretion events and VHE phenomena. 
Simultaneous data in the X-ray and soft gamma-ray ranges will 
be essential to constrain the SED and the whole picture (see 
section 4). In particular, having an important time overlap be- 
tween Fermi, X-ray satellites and CTA operations would be ex- 
tremely rewarding. This should be considered when planning 
future high-energy satellites. 

The observational study of the compact and fast-varying 
AGN and possible cataclysmic events related to SMBH physics 
will clearly benefit from an efficient alert system within a global 
network. Operating CTA under alerts will be an important strat- 
egy to catch specific AGN behaviors. This procedure is already 
followed with success by present IACT Alerts from optical 
monitoring, X-rays and gamma-rays have led to several discov- 
eries and should be pursued and extended, especially to wide- 
field instruments like HAWC and LHAASO. Thanks to its long 
lifetime, CTA will also allow to search for remarkable electro- 
magnetic counterparts of sources of non-photonic messengers 
detected by large experiments as ultra high energy cosmic rays 
with the PAO and the Telescope Array, neutrinos with IceCube 
and Km3Net and gravitational waves from SMBH with LISA 
H97l LISHl [J99l EQOl E0H . One can anticipate a great scientific 
return both in astrophysics and in fundamental physics from 
any discovery of this kind. The common occurrence of VHE 
gamma-rays and high energy cosmic rays and neutrinos is al- 
ready studied in the literature (see section 4.2), while the possi- 
ble link between VHE gamma-rays and gravitational waves still 
deserves further analysis. Non-thermal electromagnetic coun- 
terparts are expected from black hole merger (see [200| and 
references therein) and any coalescence event detected by LISA 
would be a compulsory ToO for CTA. Predictions based on the 
statistics of binary massive black holes estimate a possible rate 
of one major coalescence event every 5 years and hundreds of 
weaker events per year. However a better knowledge of the 
electromagnetic signal expected from such phenomena is still 
needed before quantifying the probability of serendipitous de- 
tection of coalescence events by CTA prior to LISA operations. 

6. Probing the intergalactic medium and backgrounds 

6.1. AGN as beacons of gamma-rays 

Beyond the studies of intrinsic AGN properties, monitoring 
of a large sample of AGN at different redshifts will provide con- 
straints on cosmic background radiation, some cosmological 
matters, and fundamental physics of space-time. Indeed, many 
observing programs related to AGN will have to be thought of 
as large key projects with multiple aims, pursuing different sci- 
entific goals although requiring almost exactly the same type of 
data. 



Photons from AGN with sufficiently high energies can pair- 
produce with photons from the infrared, optical and/or ultra- 
violet EBL in processes of the type yxev + Tir/o/uv — > <? + e~ 
(27J ESI- The IR/optical/UV EBL results from direct and re- 
processed stellar (and AGN) emission at all redshifts. In the 
early Universe, decays of exotic particles or Population III stars 
may have made substantial contributions to the EBL. The fact 
that the EBL contains calorimetric information about the en- 
tire history of IR-UV emission in the Universe makes it an ex- 
tremely interesting phenomenon for detailed studies. Although 
extragalactic pair-production processes reduce the visibility of 
high-redshift objects, studying the resultant attenuation features 
as a function of photon energy and redshift opens up exciting 
possibilities to determine the evolving intensity and spectrum 
of the IR-UV EBL, and in turn the cosmic history of star and 
galaxy formation in the Universe B202L This alternative ap- 
proach is complementary to direct measurements of the EBL, 
for which a major uncertainty is the subtraction of intense fore- 
ground components due to the zodiacal light. Measuring the 
EBL with gamma-ray absorption can circumvent this problem, 
and moreover can probe its time evolution through observations 
of sources at different redshifts. High quality spectra up to 100 
TeV and analysis of spectral features such as cut-off and absorp- 
tion features at extreme VHE will be mandatory in this regard. 
Identifying the stationary spectral pattern from variable features 
should help to disentangle EBL effects from intrinsic phenom- 
ena in the compact sources (see article on EBL, this issue). 



Monitoring of flaring AGN should also allow to check the 
Lorentz invariance and to put various limits on models of quan- 
tum gravity as VHE gamma-rays with very short wavelengths 
become sensitive to the microstructure of space-time. Tight 
constraints on relativity violations could be obtained while 
looking for a frequency-dependent velocity of gamma-rays due 
to vacuum dispersion (see the paper on Fundamental Physics, 
this issue). If present, such effects induce a possibly detectable 
time delay between arrival times of VHE flares seen at differ- 
ent energies 112051 12041 12051 [2061 12071 . They also modify the 
threshold condition for interaction of VHE particles with the 
EBL. Recent analysis of current IACT data did not allow any 
firm conclusion on this issue [138 141, 206, 207]. Such inves- 
tigation would deeply benefit from the advent of VHE detectors 
with higher level capabilities. However any time lags possibly 
detected by temporal analysis of light curves will remain diffi- 
cult to interpret, since it requires a good understanding of the 
source itself to disentangle any signature of vacuum dispersion 
along the line of sight from intrinsic time delays due to par- 
ticle acceleration, emission and radiation transfer in the AGN 
II 1391 1 1401 . Long term monitoring of a large sample of variable 
targets at various redshifts and collection of high quality light 
curves and spectra during ToO bright flares should be sched- 
uled probably as a key-programme, to optimize the probability 
of such time delay detection. Such type of multipurpose pro- 
gramme will be very rewarding in term of scientific return ver- 
sus observing time since the same set of data can address many 
different issues (AGN physics, EBL, fundamental physics ...). 
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6.2. Pair halos, pair echos and intergalactic magnetic fields 

High-energy gamma-rays from AGN offer a unique poten- 
tial to probe weak intergalactic magnetic fields (IGMF), which 
represent a long standing question of astrophysics, astroparticle 
physics and cosmology. There are various ways in which weak 
ubiquitous magnetic fields may have been generated in the early 
Universe. Such cosmological magnetic fields are potentially 
important for the description of the primordial universe, and 
could be the ultimate origin of the magnetic fields seen today 
in galaxies and clusters of galaxies by serving as the seed fields 
for subsequent amplification by galactic dynamo mechanisms 
1 208 2091 12101 . In some regions such as intergalactic voids, 
they may have survived to the present day as IGMF without be- 
ing affected by later magnetization from astrophysical sources 
B211I . and therefore may provide us with valuable information 
about physical processes in the early Universe. So far, vari- 
ous mechanisms have been proposed for the generation of such 
cosmological magnetic fields, including different types of cos- 
mological phase transitions [209, 2101 or processes related to 
cosmic reionization 1212112131 . with predicted field amplitudes 
in the range B ~ 10~ 25 - 10~ 15 G. While they may suffice as 
seeds for galactic dynamos, such tiny magnetic fields are ex- 
tremely difficult to put in evidence through conventional meth- 
ods such as Faraday rotation measurements or by their effect on 
cosmic microwave background anisotropies. 

A very powerful probe of such weak IGMF may be offered 
by the secondary GeV-TeV components accompanying the pri- 
mary TeV emission of blazars, which result from IC emission 
by e~e + pairs produced via intergalactic yy interactions among 
primary TeV and EBL photons. Depending on the IGMF value, 
such secondary components may be observable either as "pair 
echos" that arrive with a time delay relative to the primary emis- 
sion I214, l215ll216l l2261. or as extended emission with a spatial 
extension around the primary source M217II2181I219II . The prop- 
erties of the extended emission depend on the IGMF strength. 
Strong enough IGMF leads to full isotropization of the cascade 
emission and formation of a physical "pair halo" around the 
primary gamma-ray source B217r | while weak magnetic field 
leads to appearance of an extended emission with IGMF depen- 
dent size, which emission benefits from Doppler boosting. 

There have been no clear observational indications for either 
pair echo or halo emission so far, and deriving lower bounds 
to IGMF strengths from upper limits to such components ob- 
tained by GeV-TeV observations of selected blazars appears 
very promising and still a matter of debate [ 220 221 , 222 , 223 
224, 225 226 1 . Existing limits rely on a combination of data 
obtained with different types of telescopes in the GeV and TeV 
bands. CTA will have the great advantage of being able to 
cover all the relevant energy range, from 10 GeV up to 10 TeV, 
within one observational facility. Improving the time coverage 
to better determine the GeV and TeV light-curves with Fermi 
and hopefully future GeV gamma-ray satellites, together with a 



global network of IACT including the two sites of CTA, would 
significantly benefit to this issue. 

The positive detection of the pair echo or extended emis- 
sion components would provide highly valuable information on 
the IGMF, a goal that is expected to be achievable with CTA. 
Pair echos could be detected in the case of very tiny IGMF 
(< 10~ 16 G), and extended emission for higher IGMF values. 
CTA will conduct detailed measurements of spectral variability 
with higher sensitivity and time resolution over a wider energy 
band compared to current gamma-ray telescopes, potentially al- 
lowing us to disentangle and positively identify the echo com- 
ponent from the primary emission, and thereby probing very 
weak IGMF in the range B < 10~ 16 G, totally out of reach with 
other means of investigation. 

A fruitful synergy with wide-field instruments such as 
HAWC or LHAASO can be expected in the CTA era, since reg- 
ular, long-term coverage of the multi-TeV emission of blazars 
by such facilities should greatly improve our knowledge of their 
TeV activity on timescales of years to decades, giving us a much 
better handle on the light curve of the primary emission as in- 
put for more reliable predictions of the pair echo and extended 
emission properties. The prediction of the level of the cascade 
emission will be greatly improved also with the detailed mea- 
surements of the blazar spectra in the 10 TeV band with CTA 
itself, which will give the details of absorption of the primary 
gamma-rays in the intergalactic medium. 

If IGMF strengths are in a stronger range, B ~ 10~ 16 - 
10 12 G, the spatially-extended emission may be detectable and 
resolvable by CTA by virtue of its high sensitivity and angular 



resolution (see Fig. 19 1, in which case one would probe not only 
IGMF but also the lifetime of TeV activity for the source. High 
signal-to-noise detections of pair halos around nearby blazars 
(z ~ 0.1 - 0.2) would allow us to study the energy -dependent 
radial profile and thereby deduce unique information on the 
IGMF and on the EBL (228). 

If the IGMF is B > 10~ 12 G and strong enough to effectively 
isotropize the secondary pairs, the detection of isotropic pair 
halos from distant AGNs should be feasible [217, 230 1. Figure 



20 shows the expected sizes of such isotropic pair halos as a 



For an isotropic pair halo to form, there must be a relatively high mag- 
netic field (B > 10 G) within ~ 10 Mpc of the source to facilitate isotropic 
emission. 



function of redshift, generated by 10 TeV primary photons giv- 
ing rise to ~100 GeV secondary photons, and using the EBL 
model of Franceschini et al. 2008 fiTTl . Even though the prop- 
erties of such isotropic halos no longer depend explicitly on the 
IGMF, observation of their isotropic nature would still provide 
indicative lower bounds on its value, allow to check the self- 
consistency of our EBL description, and also give important 
information on the high-energy spectrum and beaming angle of 
the primary emission. 

Although the halos are very large for redshifts z < 0.5, a halo 
with a radius < 1° at z > 0.5 would fit into the field of view of 
CTA telescopes and would allow for reliable background sub- 
traction. If one assumes that the source has a constant vF v -flux 
from 100 GeV to 10 TeV, a simple estimate shows that the halo 
fluxes are below a couple of percent compared to the primary 
100 GeV emission. The halo fluxes would be much larger if the 
photon index at the source was much harder than r = 2. This 
is to be compared to the pair halo detection limits which should 
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Figure 19: The arrival directions of the primary and secondary gamma-rays 
(circles) from a source at a distance D = 120 Mpc. The IGMF strength is 
10~ 14 G (upper panel) and 10~ 15 G (lower panel). The sizes of the circles rep- 
resenting each photon are proportional to the photon energies. The AGN intrin- 
sic gamma-ray spectrum is typically described as dN/dE ~ E 2 e.xp(-E/E cut ) 
with the cutoff energy E cul ~ 300 TeV. The blue dashed and red solid circles 
are radii of 1.5° and 2.5°, respectively, and would perfectly fit into the FoV of 
CTA. See 1229] for details. 



be actually reachable with CTA. 

The expected flux sensitivity for pair halo detection with 
CTA has been calculated for four separate approaches. In all 
cases a differential angular distribution of a pair halo at z - 



0.129 and E y > 100 GeV, taken from Fig. 6 of I23TI was used 
as our model halo profile, with a functional form: dN / df? 2 oc 
9 5 ^. The differential energy flux sensitivity for the four differ- 
ent approaches is shown in Fig 21 Methods A and B rely on 



a 5 cr excess above the expected background, calculated using 
Equation 17 in [61 1. Method A searches for some overall ex- 
tended emission with a halo profile for 6 < 0.32°. Method B 
probes a region in which a point-like central source would no 
longer be dominant (0.11° < 9 < 0.32°) and therefore provides 
the most basic method of establishing the detection of extended 
emission. For methods C and D the "goodness of fit" for a halo 
profile (convolved with the CTA PSF) fitted to simulated CTA 
data is used to determine the expected flux sensitivity. Method 
C compares a halo fit to the null hypothesis of background flue- 





Figure 20: Radius in Mpc (upper panel) and angular radius in degree (lower 
panel) at which the y-ray emissivity of isotropic pair halos drops by a factor of 
1/e as functions of redshift. The calculation assumes 10 TeV primary photons 
leading to 100 GeV secondary photons in the observer frame, an IGMF strength 
strong enough to isotropize the secondary pairs, and the EBL model of 1411 . 



tuations only, and the optimal 6 limit for the fitting was found 
to be « 0.2°. Method D however, tests how well a point-like 
source and a pair halo can be distinguished. This is done by 
assessing the difference between the likelihood obtained for a 
PSF hypothesis and that obtained for a pair halo function, with 
the limiting flux defined as that at which the point-source hy- 
pothesis can be rejected at the 5 cr level. The rather high flux 



sensitivities for Methods B and D in Fig 21 indicate the diffi- 
culty of identifying halo-like emission at this confidence level 
(CL). If a 95 % CL is taken instead of 5 cr, then the sensitivity 
of Method D is comparable to Method A and C. Therefore a po- 
tentially extended source detected at the sensitivity limit with a 
significance of 5 cr could be distinguished from a point source 



with about 95% confidence. Note that in Fig 21 the flux quoted 
is that within a 1° region of the source. The lower panel of 



Fig 21 shows the pair halo sensitivity using Method A derived 
for a range of different CTA candidate configurations. Config. 
I appears to be the most suitable setup for pair halo studies. 

To summarize, any detection of pair halos, extended emis- 
sion, or echos with CTA would be a great step forwards in the 
characterization of cosmic magnetic fields and backgrounds. It 
could also open new paths to cosmologists and astroparticle 
physicists, if constraints on the primordial magnetic field could 
be deduced from the IGMF properties [ 209 1 . 

6.3. Observing the diffuse VHE background with CTA 

Extragalactic diffuse background radiation has been ob- 
served from radio to GeV energies. Evidence for diffuse extra- 
galactic gamma-ray background (EGRB) has been reported first 
by the SAS 2 satellite £521 . and more firmly by the EGRET de- 
tector [ 233 234 1 . The separation of an extragalactic component 
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Figure 21: Upper panel : A comparison of the sensitivity of CTA configura- 
tion I to pair halo emission for four different analysis methods. These methods 
are described in the main text. Lower panel : Flux estimates on the expected 
pair halo emission with CTA, for various CTA array configurations (50 hours 
observing time for each field, 20 degrees zenith angle observations), using anal- 
ysis Method A. A differential angular distribution of a pair halo at z = 0. 1 29 and 
E y > 100 GeV, taken from Fig. 6 in [23 1] was used as our theoretical model 
(dashed line). 



from foreground and galactic emission is generally challenging 
and model-dependent, and the detection of EGRB was a mat- 
ter of debate |235|. The origin of EGRB is still unsolved, but 
there must be some contribution from unresolved or unidenti- 
fied astrophysical sources, especially blazars, active nuclei and 
galaxies (see Fig 22 1, More exotic processes like <? + e~ pair ha- 



los or even self-annihilating dark matter could also contribute. 
Recently Fermi reported their new results about EGRB 

1244, 2451 



The reported flux is consistent with the EGRET 
data at 100 MeV, but the Fermi spectrum is softer than EGRET 
and the Fermi flux is lower than EGRET by a factor of 2-3 
at 10 GeV. The Fermi team suggested that unresolved blazars 



can explain at most ~23% of the Fermi EGRB flux at > 100 
MeV, indicating a significant contribution from other popula- 
tions. Indeed, the Fermi EGRB spectrum is well described by 
a power-law with photon index of y — 2.4, which is in good 
agreement with the prediction based on theoretical models of 
EGRB from astrophysical sources like blazars, radio galax- 
ies, and star- forming galaxies ll55l 12461 12471 . Against the to- 
tal EGRB flux (unresolved diffuse plus bright resolved sources, 
~ 1.4 x 10~ 5 cm^s -1 at > 100 MeV), the predictions by these 
models are 44% by blazars, 18% by radio galaxies, and 7% by 
star-forming galaxies. Therefore, >70% of the total EGRB flux 
can be explained by the known population of gamma-ray emit- 
ters. The remaining 30% may come from new, still undetected 
populations, or from systematic uncertainties in the EGRB 
measurement, e.g., foreground subtraction. Fermi EGRB mea- 
surement below 100 GeV has already set an upper-limit on the 
EGRB itself above 100 GeV [248 1. The limit is conservative for 
the cascade emission from EGRB interacting with the cosmic 
microwave-to-optical background radiation not to exceed the 
current EGRB measurement. The cascade component fits the 
measured VHE EGRB spectrum. However, taking into account 
known components, the upper limit can be approximated as 
E 2 dN/dE < 4.5 x 10- 5 (£/100GeVr ' 7 MeV/cm 2 /s/sr [248] 
which is below the Fermi EGRB measurement above 100 GeV. 
This suggests the existence of new physics in this energy range 
and shows the need to study the VHE EGRB for a better under- 
standing of the VHE gamma-ray sky. 

Detection of the VHE gamma-ray background would provide 
important information to solve the EGRB problem. Especially, 
we expect a turnover above 100 GeV due to intergalactic ab- 
sorption by the EBL, and the break energy depends on the red- 
shift distribution of the sources contributing to EGRB. More- 
over, HBL population like 1ES 0229+200 are not detected by 
Fermi but by current IACT However, it is challenging to de- 
termine the diffuse background level by an observation using a 
ground-based Cherenkov telescope, because the photon events 
must be discriminated from the much more frequent cosmic- 
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ray electron events. In the measurement of cosmic -ray flux, the 
H.E.S.S. team found that the distribution of the depth of shower 
maximum (X m - dx ) indicates less than 10% of photon contribu- 
tion in the electron flux, though they also mentioned that ~50% 
photon contribution cannot be excluded when the systematic 
uncertainty in X m . dx determination is conservatively taken into 
account. The 10% level of electron flux is indicated in Fig 22 
The expected EGRB signal in the TeV range is about 1% of 
the electron flux. Detecting it with CTA will be challenging, 
but if it is achieved, we will get an important insight into the 
extragalactic VHE populations. 

7. Conclusion 

The discovery of AGN at extreme energies was a real break- 
through performed by IACT of the 90s decade. The current 
IACT increase the VHE extragalactic sample, and provide ba- 
sic spectral and variability characteristics. They start constrain- 
ing the physical processes at work and to catch a glimpse of 
the extragalactic VHE populations. Together with EGRET and 
later on Fermi, they show that the emitted gamma-ray power 
can represent a significant part of the non-thermal activity en- 
ergy budget, and they raise a lot of intriguing questions. 

CTA will have the unique potential to open new science di- 
rections beyond the conventional AGN paradigm, and to shed 
new light on AGN and SMBH physics, populations, and evo- 
lution. Especially the jump achieved in sensitivity and spec- 
tral resolution, and the large spectral coverage from a few tens 
of GeV to several tens of TeV will be mandatory in this re- 
gard. Improvement of the angular resolution and of astromet- 
ric accuracy will facilitate source identification and may lead 
to the discovery of diffuse or multiple VHE components asso- 
ciated to AGN. The numerous modes of observation foreseen 
for the CTA array are perfectly appropriate to AGN observing 
programs and should definitely be implemented. Needless to 
say that AGN studies will strongly benefit from having access 
to the whole sky. Developing the two sites of the CTA array, 
one in the southern and one in the northern hemisphere, will be 
mandatory in this regard. Having them at different longitudes 
would allow to improve the time coverage for variable sources 
reachable from both sites. As far as the current CTA candidate 
arrays are concerned, emphasis on LSTs could provide a cru- 
cial boost in the effective area in the low-energy region where 
high-redshift sources are expected to be most prominent. So- 
lutions without LSTs would have a negative impact on AGN 
population studies and every effort should be made to smoothly 
connect the low-energy regime of CTA (< 30 GeV) with Fermi. 
Conversely, the VHE fluxes of most of the bright Fermi blazars 
are largely above the CTA sensitivity limit and interesting high- 
statistics spectral information on these blazars will be available. 
In particular, observing over a large spectral range up to sev- 
eral tens of TeV with a good spectral resolution will make it 
possible to find out whether the observed cut-offs in the blazar 
spectra are intrinsic to the source or are induced by the effect of 
EBL absorption, and to analyze the maximal energies of parti- 
cles responsible for the VHE emission and extreme acceleration 
processes. 



Hundreds of VHE AGN will be detected over the years of 
CTA operations. It will be the role of CTA to produce re- 
fined light curves to probe the shortest time scales, high quality 
spectra to constrain acceleration, emission and absorption ef- 
fects, and maps to locate the VHE emission zones of nearby 
radiogalaxies. AGN monitoring, serendipitous discoveries and 
deep surveys will gather meaningful samples of VHE sources to 
contribute to the elaboration of a convincing unifying scheme 
for AGN. AGN physics appear as a science driver for CTA, 
exploring present mysteries such as mechanisms to accelerate 
VHE particles in extremely compact zones. Generally speak- 
ing, VHE instruments offer a completely new and independent 
way to explore AGN and their most efficient and energetic re- 
gions. They should provide a decisive key for the understand- 
ing of the origin of several non-thermal phenomena at work in 
such sources, in particular the missing link between accretion 
physics, SMBH magnetospheres and jet formation and propa- 
gation, which have been studied for several decades at lower 
energies and still remain fundamental open questions of con- 
temporaneous astrophysics. 

As a caveat, it should be noted that all quantitative estimates 
presented in this article are preliminary, as further refinements 
to the CTA array configurations and improvements to the anal- 
ysis tools are being worked on. Uncertainties in the spectral pa- 
rameters and the EBL model could introduce biases. However, 
it is clear that we can look forward to a densely populated sky 
in the TeV range with CTA. This will open new opportunities 
to extragalactic studies, including dedicated multiwavelength 
campaigns and studies of very rapid flares. Furthermore, CTA 
will matter for research on blazar evolution, with the detection 
of high-redshift blazars. 

CTA will be a crucial instrument for multi-messenger as- 
tronomy. Indeed, hadronic scenarios for TeV emission pre- 
dict the simultaneous production of neutrinos and cosmic rays. 
Bright objects detected by CTA will be prime candidates for 
non-photonic sources potentially reachable by large cosmic ray 
and neutrinos detectors such as the PAO and the KM3Net. Con- 
versely, any future detection of cosmic non-photonic sources of 
neutrinos, cosmic rays or even gravitational waves will benefit 
from follow-up observations with CTA to identify their nature, 
further explore their properties and constrain the parameters. 

While observing AGN with CTA, many other questions less 
directly related to AGN topics, can be simultaneously explored, 
using AGN as bright remote beacons to analyse various proper- 
ties of the intervening media, such as the extragalactic back- 
ground light, and the incredible perspective to measure ex- 
tremely tiny values of the intergalactic magnetic field, or even 
the possible quantum properties of space-time. Unexpectedly, 
CTA turns out to be an ideal instrument to probe the nature of 
the seed IGMF and could solve the problem of the origin of 
magnetic fields in galaxies and galaxy clusters which is one of 
the long-standing unsolved problems of astrophysics and cos- 
mology. 

An important synergy between the analysis of the future 
VHE data and further developments of theoretical related fields 
is anticipated and will be mandatory to handle the flow of new 
information expected from the CTA observatory. 
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Appendix A. Appendix on extrapolation of Fermi spectra 

For the conservative estimate, based only on AGN with con- 
firmed redshifts in the Bzcat or Veron (13th edition) catalog and 
without any flags in the 2FGL catalog, 508 sources from the 
cross-correlated 2FGL catalog were selected for further pro- 
cessing, including 288 flat-spectrum radio quasars (FSRQs), 
156 BL Lacs, 7 radio galaxies (RGs), 5 Seyfert galaxies, 2 star- 
burst galaxies (SBGs) and 50 AGN of other classes. In the sec- 
ond study, for which additional sources of redshift estimates 
where used and the 2FGL analysis flags where ignored, a sub- 
set of 561 Fermi sources was selected, including 340 FSRQs, 
171 BL Lacs, 10 RGs, 6 Seyfert galaxies, 3 SBGs and 31 AGN 
of other classes. 



To extrapolate the spectra to the highest energies, we use 
the integral flux from Fermi measurements between 1 and 
100 GeV ("F1000") or the differential flux at the pivot energy 
("FluxJDensity"), together with the spectral index (F) given in 
the 2FGL catalog. For each individual source, we adopt the 
corresponding power-law or LogParabola parametrization as 
described in |59l . For very hard sources following a power- 
law with F < 2 in the Fermi range, a straight extrapolation 
could create runaway integrations, therefore we apply an arti- 
ficial break in the power law with a photon index of Y2 — 2.5 
above 100 GeV to soften such spectra. The latter is in broad 
agreement with observed AGN features. 

Using the extrapolated and absorbed AGN spectra, we in- 
tegrate the flux of particles per energy bin weighted with the 
effective areas for the different candidate arrays. The resulting 
expression is subsequently multiplied by the observation time, 
yielding the total number of detected photons. For each source, 
the significance was calculated using Equation 17 in [61] as- 
suming N on (on region) to be the number of source photons 
plus the number of photons derived from the simulated back- 
ground rates, and N s fixed at the background rate (off-region). 
The variable a is given by the ratio of the sizes of the two re- 
gions times the ratio of the exposure times and the respective 
acceptances. 

For the first study, a is assumed to be 0.1. The detection of 
a source requires a significance of the signal over noise ratio 
of at least 5<x, a minimum of 7 excess events and a signal ex- 
ceeding 3% of the background. H For the second study, an a 
of 0.2 is used and a signal over 5% of the background is re- 
quired. To maximise the number of detections in this study the 
energy threshold for the detection is varied such as to optimize 
the signal from each source. 

Appendix B. Appendix on variability statistical analysis 

The unbinned light-curves were generated using the algo- 
rithm of [ 181 1, as a sequence of time-tagged events, which can 
then be directly inspected for variability with no recourse to 
binning. To analyze the variability properties of the unbinned 
datasets in a non-parametric fashion, we adopted a Bayesian al- 
gorithm called "Bayesian blocks". The method uses a dynamic 
algorithm to analyze the entire light-curve and find the optimal 
way of dividing it into a series of contiguous blocks, k, each 
consistently describable as an homogeneous Poisson process of 
rates A^. The change-point between each block marks the po- 
sition where a significant change in the rate of arrival of events 
occurred, the entire event-sequence being thus treated as an in- 
homogeneous Poisson process made up from the juxtaposition 
of constant-rate processes. The approach is therefore very con- 
venient as a model-independent analysis of the time-series, in 
the sense that the partitioning of the light-curve is based only on 
the (trivial) hypothesis that the event sequence follows a Pois- 
son distribution. At every time in which the hypothesis of con- 
stant rate for the distribution is no longer valid for a subsection 

The CTAmacrosv5 and CTAmacrosv6 tools for the estimation of the in- 
strumental response were used for this study. 
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of the events, a new partition (or block) is created, within which 
the approximation of homogeneity is again valid. In this way 
the partitioning naturally selects the flares in the time series and 
the variability timescale of a given feature can be directly iden- 
tified with the size of the corresponding block. If one wants to 
quote the variability in function of the more common quantity 
of the "doubling time" Td of the flux, a simple correspondence 
can be made, by which the size of the block Aty can be taken 
as approximately 2<x of a Gaussian fit to the flare events. In the 
case of generalized Gaussian profiles, such as presented in the 
text, the correspondence can be made in terms of Aty ~ cr r + <x d , 
where the different rise and decay times of the profile are taken 
into account. For a Gaussian profile, the conversion factor is: 
Af k = 2/ln(2)cr * 3r d . 

As with any Bayesian method, an a priori hypothesis on the 
time series model must be tested against the data. In the case of 
the Bayesian blocks, the prior is on the number of blocks which 
will compose the partition model of the event sequence. With- 
out entering into much detail, the algorithm uses a geometric 
prior which penalizes for excessive division of the light-curve 
into too many parts, as it is likely that the number of blocks 
(flares), k, is k <c N, the number of events in the sequence. 
For each new block added in the light-curve, a penalizing term 
ln{y) is added to the log-likelihood of the partition model. The 
value of the prior y here is a free parameter to be determined by 
simulation, but a strength of the model is that for a large range 
of values of y, the model finds a stable, optimal solution, as can 



be seen in Fig B.23 For the analysis shown in the text, a prior 
of ln(y) = 6 was used, corresponding to a probability odds for 
block division of ~ 5 x 10 2 . 



The results presented in Fig B.23 were derived by simulating 
two groups of 1000 light-curves of about 10 3 and 10 4 events 
each, for H.E.S.S. and the CTA extrapolation. The two curves 
shown present the mean block size Aty and the variance of 
the histogrammed distributions for both sets of light-curves: 
H.E.S.S. -like (upper-curve) and CTA-like (lower-curve). The 
plateau part of the curve indicates the regions where the al- 
gorithm gives stable solutions to the partitioning of the light- 
curve, chosen for the analysis, the gray-bands on the figure 
indicate the \<j levels of the histogrammed distributions of the 
block sizes. 
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Figure B.23: Analysis of the stability of the prior to the partition model, In(y). 
The plot shows the distribution of the sizes of the recovered blocks versus 
the prior values. As can be seen, for too low priors, the light-curve is over- 
divided, as there is little penalization for the modeling of the series with too 
many blocks. Conversely, for prior values greater than \n(gamma) ~ 8, the 
number of blocks starts to slowly decrease (sizes of blocks increase), indicat- 
ing that the penalization introduced by the prior is now strongly influencing the 
final partition. Between 3-6 the prior shows an optimal, stable behaviour. 
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